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Summary 


In  the  fourth  reporting  period  we  have  extended  and  concluded  the  theoretical  analysis  of  noise 
propagation  in  OPCPAs,  and  gained  useful  information  for  the  optimization  of  OPCPA  design 
parameters.  From  the  experimental  side,  we  have  optimized  two  set-ups  for  coherent  pulse  synthesis 
and  CEP-stable  pulse  generation. 

The  main  results  obtained  are  summarized  in  the  following: 

Simulations  of  noise  evolution  in  OPCPAs 

The  quantum-mechanically  consistent  numerical  model  detailed  in  the  previous  reports  allowed 
us  to  conclude  the  investigation  of  the  dynamics  of  superfluorescence  growth  in  a  realistic  high-gain 
OPCPA.  Since  the  model  can  simulate  also  the  saturation  stage  of  amplification,  our  investigation  for 
the  first  time  captured  all  dynamics  of  a  quantum-noise-contaminated  OPCPA,  addressing  the 
important  practical  issues  of  signal  energy  stability  and  pulse  contrast.  The  different  evolution 
dynamics  of  these  quantities  were  analyzed  throughout  the  amplification  process  in  three  operating 
conditions,  characterized  by  different  chirps  of  the  input  seed.  We  find  maximization  of  the 
efficiency-bandwidth  product  is  correlated  to  the  noise  contamination;  in  addition,  while  amplifier 
saturation  improves  the  signal’s  shot-to-shot  energy  stability,  it  does  not  necessarily  improve  the  pulse 
contrast.  Knowledge  of  these  dynamics  increases  our  fundamental  understanding  of  quantum  noise  in 
parametric  amplification;  it  also  provides  important  insight  for  the  optimization  of  OPCPA  systems 
applied  to  the  study  of  strong-field  laser  physics. 

Synthesis  from  two-stage  parametric  amplification 

In  order  to  generate  the  gap-free  phase-stable  ultrabroadband  WLC  required  by  this  task,  we  have 
first  developed  a  two-stage  IR-OPA.  This  setup  removes  the  spectral  gap  at  800  nm  typically  arising 
after  spectral  broadening  the  fundamental  frequency  of  a  Ti:sapphire  laser.  The  two-stage  IR-OPA 
provides  signal/idler  pulses  at  wavelengths  longer  than  1.3  pm,  which  are  then  broadened  in  a  YAG 
or  sapphire  plate.  The  WLC  is  split  in  two  and  used  to  seed  simultaneously  the  visible  NOPA  and  the 
degenerate  OP  A,  which  are  then  separately  compressed  using  chirped  mirrors  to  nearly  TL  duration. 
The  two  pulses  are  combined  using  a  broadband  beam  splitter  and  their  relative  CEP  and  delay  is 
monitored  by  spectral  interferometry.  Preliminary  results  demonstrated  that  combination  of  the  OPA 
pulses  allows  obtaining  a  broadband  pulse  extending  from  550  to  1000  nm,  corresponding  to  a 
transform- limited  duration  of  4  fs. 

Sub-cycle  pulse  synthesis  from  two  OPCPAs 

We  have  realized  a  waveform  synthesis  scheme  which  is  able  to  combine  high-energy,  few-cycle 
optical  pulses  from  multi-color  OPCPAs.  The  system  generates  a  sub-cycle  waveform  with  a 
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spectrum  spanning  close  to  two  octaves  and  15-pJ  pulse  energy.  This  source  can  be  applied  to  the 
direct  generation  of  isolated  soft-x-ray  pulse  by  HHG,  eliminating  the  need  for  gating  techniques  or 
spectral  filtering.  The  system  is  capable  of  stabilizing  and  controlling  all  independent  parameters  that 
define  the  synthesized  electric-field  waveform,  such  as  the  absolute  phase  and  the  relative  delay  of  the 
two  pulses. 

Introduction 

The  aim  of  the  project  is  to  develop  two  (or  more)  ultra-broadband  CEP-stable  Optical  Parametric 
Amplifiers  (OP As)  and  coherently  combine  their  outputs  in  order  to  synthesize  high  energy 
(sub)single-cycle  pulses  with  controlled  electric  field  profile. 

OP  As  are  devices  which  exploit  second  order  optical  nonlinearity  in  order  to  efficiently  transfer 
the  energy  of  a  fixed  frequency  “pump”  pulse  to  a  broadly  tunable  weak  “signal”  pulse.  To  achieve 
efficient  energy  transfer,  the  so-called  “phase  matching”  condition  between  the  interacting  waves 
must  be  satisfied.  Under  appropriate  conditions,  phase  matching  can  be  achieved  over  a  broad 
frequency  range,  thus  turning  OP  As  into  ultra-broadband  amplifiers  and  allowing  the  generation  of 
tunable  few-optical-cycle  light  pulses. 

Broad  gain  bandwidth  in  an  OPA  is  achieved  when  the  group  velocities  of  signal  and  idler  are 
matched  [1];  this  condition  is  satisfied  either  in  the  case  of  type  I  phase  matching  at  degeneracy,  or  in 
the  non-collinear  OPA  (NOP A),  in  which  the  idler  group  velocity  is  projected  along  the  signal 
propagation  direction.  Using  these  concepts,  a  variety  of  broadband  OPA  schemes,  pumped  by  either 
the  fundamental  frequency  (FF)  or  by  the  second  harmonic  (SH)  of  Ti: Sapphire  and  seeded  by  white 
light  continuum  (WFC),  have  been  demonstrated.  The  SH-pumped  NOP  A  in  the  visible  [2-4],  the  FF- 
pumped  NOP  A  in  the  near-IR  [5,  6]  and  the  FF-pumped  near-IR  degenerate  OPA  [7]  allow  to  cover 
nearly  continuously  the  wavelength  range  from  500  to  2  pm. 

The  goals  of  this  project  are:  (i)  to  develop  two  broadband  OP  As  for  the  amplification  of  CEP- 
stable  seed  light  in  a  broad  wavelength  range,  from  500  nm  to  2  m;  (ii)  to  coherently  combine  these 
OP  As  in  order  to  synthesize  a  single-cycle  optical  pulses  with  CEP  control,  i.  e. ,  with  a  precisely 
defined  electric  field  waveform  [8].  Such  CEP-stable  single-cycle  pulses  will  be  used  to  drive  strong- 
field  processes  such  as  high-harmonic  generation  (HHG),  which  has  recently  led  to  efficient 
generation  of  isolated  attosecond  pulses,  with  duration  as  short  as  80-as  in  the  XUV  range  [9]. 

In  order  to  fulfill  the  aims  of  the  project,  the  following  two  tasks  have  been  addressed  in  the 
fourth  period  of  the  project: 
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Task  1  -  Generation  of  few-optical-cycle  CEP  stable  pulses  from  OPAs 


Simulations  of  noise  evolution  in  OPCPAs 

Our  pulse  synthesis  experiments  require  the  coherent  combination  of  pulses  generated  by 
parametric  amplification  (OPAs  and  OPCPAs).  The  unusual  phenomena  associated  with  quantum 
noise  contamination  in  an  OPCPA,  observed  but  heretofore  not  well  understood,  are  the  concern  of 
this  part  of  the  project.  While  OPCPAs  are  immune  from  amplified  spontaneous  emission,  they  are  on 
the  other  hand  prone  to  another  parasitic  effect  with  evolution  dynamics  that  are  not  well  understood 
in  this  geometry,  namely  parametric  superfluorescence  (PSF),  i.e.,  parametric  amplification  of  the 
vacuum  or  quantum  noise  due  to  two-photon  spontaneous  emission  from  a  virtual  level  excited  by  the 
pump  field  [10].  Previous  experimental  results  have  shown  that  the  presence  of  PSF  results  in  an 
amplified  signal  field  with  two  apparent  macroscopic  components:  a  coherent  pulse  with  well  defined 
temporal  chirp  that  matches  that  of  the  seed  pulse,  and  which  is  dechirped  in  the  compression  stage  to 
generate  a  transform-limited  pulse,  and  an  incoherent  pedestal  with  stochastic  phase  statistics  similar 
to  that  of  spontaneous  parametric  generation,  which  cannot  be  recompressed.  Henceforth,  we  refer  to 
these  phenomenological  components  observed  at  the  output  of  an  OPCPA  as  the  coherent  pulse  and 
incoherent  pedestal ,  respectively.  Here  we  will  report  on  the  numerical  study  of  the  incoherent 
pedestal;  the  aim  of  the  study  is  to  provide  a  guideline  for  the  project  and  realization  of  OPCPA  with 
optimized  signal-to-noise  ratios. 

Task  2  -  Coherent  pulse  synthesis 

We  have  implemented  two  schemes  for  the  synthesis  of  nearly  single  cycle  pulses:  coherent 
combination  of  two  broadband  OPAs  (performed  in  Milano)  and  of  two  high-energy  OPCPAs  in  the 
near  and  mid-IR  (performed  at  MIT).  The  first  scheme  offers  the  advantage  of  a  clean  gap-free 
spectrum,  allowing  the  generation  of  pulses  without  satellite  lobes,  but  has  currently  limited  output 
energy;  the  second  scheme,  on  the  other  hand,  has  the  advantage  of  energy  scalability  but  presents 
some  side-lobes  due  to  a  gap  in  the  synthesized  spectrum. 

Synthesis  from  two-stage  optical  parametric  amplification 

In  the  previous  report  we  introduced  three  possible  schemes  for  the  synthesis  of  few-cycle  pulses, 
based  on  the  combination  of  two  ultrabroadband  OPAs  powered  by  the  same  Ti:  Sapphire  laser.  OPAs 
are  very  versatile  and  efficient  tools  for  the  generation  of  few-cycle  optical  pulses.  In  this  phase  of  the 
project,  we  are  developing  the  first  of  the  proposed  schemes,  based  on  the  coherent  combination  of 
two  separate  OPAs  operating  in  parallel. 
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Sub-cycle  pulse  synthesis  from  two  OPCPAs 

We  demonstrate  a  new  approach,  based  on  coherent  wavelength  multiplexing  of  ultra-broadband 
OPCPAs,  for  the  generation  of  fully  controlled  high-energy  sub-cycle  optical  waveforms  with  spectra 
spanning  close  to  two  octaves.  Such  pulses  can  be  used  to  efficiently  generate  isolated  attosecond 
pulses  without  the  use  of  gating  techniques  [11].  The  system  coherently  combines  two  CEP- 
controlled,  few-cycle  pulses  obtained  from  different  OPCPAs:  1)  a  near  infrared  (NIR)-OPCPA, 
producing  25-pJ,  9-fs  pulses  centered  at  870  nm;  and  2)  a  short-wavelength  infrared  (SWIR)-OPCPA, 
producing  25-pJ,  24-fs  pulses  centered  at  2.15  pm. 

Methods ,  assumptions  and  Procedures 

Task  1  -  Generation  of  few-optical-cycle  CEP  stable  pulses  from  OPAs 
Simulations  of  noise  eyolution  in  OPCPAs 

The  quantum-mechanically  consistent  numerical  model  already  introduced  in  the  previous  reports  was 
further  optimized  and  applied  to  evaluate  the  nature  of  noise  in  OPCPAs.  Such  noise  arises  from 
parametric  superfluorescence,  i. e.  amplification  of  vacuum  fluctuations.  For  the  numerical  description 
of  this  process,  we  focus  on  an  OPCPA  seeded  by  the  initial  quantum  noise  field  and  a  chirped  signal 
field.  The  nonlinear  quantum  system  dynamics  can  be  described  by  a  quasi-probability  distribution, 
such  as  the  Wigner  distribution  (WD)  [12].  For  linear  systems  the  evolution  equation  for  the  WD  is 
equivalent  to  a  classical  Fokker-Planck  equation,  and  is  thus  also  equivalent  to  a  stochastic  process 
involving  classical  noise  sources,  resulting  in  a  semiclassical  picture  of  the  quantum  process.  This 
correspondence  has  been  exploited  in  numerical  studies  of  PSF,  OP  A,  and  optical  parametric 
oscillation  in  their  linear  regimes  [13,  14].  The  Fokker-Planck  approximation  also  holds  for  the  case 
of  weak  nonlinearities,  and  the  nonlinear  quantum  system  dynamics  can  still  be  extracted  accurately 
from  stochastic  Fangevin  equations  [15,  16],  an  approach  used  earlier  to  study  the  quantum  noise  in 
parametric  amplifiers  used  for  squeezed  light  generation  [17,  18].  These  stochastic  equations  have  a 
deterministic  component  equal  to  the  Heisenberg  equations  of  motion  for  the  field  operators  and  are 
complemented  by  relaxation  terms  and  associated  noise  terms.  For  the  case  of  a  lossless  OPA  process, 
fluctuations  stem  solely  from  the  quantum  mechanical  uncertainty  in  the  input  fields.  Knowledge  of  a 
quasi-probability  distribution  allows  computation  of  all  expectation  values  of  quantum  mechanical 
observables,  and  for  the  case  of  the  WD  and  its  associated  stochastic  process,  computed  expectation 
values  correspond  to  quantum  mechanical  expectation  values  of  symmetrically  ordered  field  operators 
[12].  In  practice,  we  simulate  the  evolution  of  noise  in  an  OPCPA  by  numerically  solving  the  coupled 
nonlinear  equations  of  parametric  amplification  in  the  spectral  domain,  accounting  for  linear 
dispersion  to  all  orders  [19,  20].  Our  1-D  plane  wave  model  includes  all  longitudinal  modes,  m,  and 
their  associated  noise.  In  the  frequency  domain,  at  any  mode  frequency,  com,  the  corresponding 
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component  of  the  initial  signal,  idler,  or  pump  electric  field  is  represented  by  a  complex  stochastic 
phasor,  Am(0)=Bm+^m.  Bm  is  the  deterministic  component  of  the  field,  and  is  set  to  0  in  the  case  of  the 
idler;  nm  is  a  zero-mean,  stochastic  phasor  representing  the  independent  fluctuations  of  the  field.  The 
initial  fields  then  follow  the  conventional  interaction  and  propagation  equations  of  the  OPCPA 
process.  Note  that  fluctuations  are  included  for  each  of  the  signal  and  idler  fields.  Real  and  imaginary 
components  of  nm  are  taken  as  uncorrelated  Gaussian  distributions  [13]  with  variance  am2  oc  <Dm. 

By  virtue  of  the  model’s  adherence  to  quantum  mechanics  even  in  the  saturation  stage  of 
amplification,  our  investigation  captures  for  the  first  time  the  saturation  dynamics  of  quantum-noise- 
contaminated  OP  A,  a  problem  of  general  interest  in  the  fields  of  quantum  and  nonlinear  optics,  both 
theoretically  and  experimentally. 

We  have  simulated  the  broadband  2-pm  OPCPA  which  has  been  developed  in  the  labs  at  MIT. 
The  system  employs  a  periodically -poled  lithium  tantalate  crystal  pumped  by  a  9-ps  Gaussian  pulse  at 
1.047  pm  and  seeded  by  a  broadband  pulse  at  2.094  pm,  for  operation  around  degeneracy.  The  seed 
spectrum  has  FWHM  bandwidth  of  69  THz,  which  well  matches  the  15-fs  phase-matching  bandwidth 
of  the  amplifier.  We  investigated  a  high-gain  amplifier  stage,  where  the  pump-to-signal  energy  ratio  is 
Ep/Es=  106.  The  seed  field  was  chirped  respectively  to  1.0-ps,  7.3-ps,  and  10.5-ps  durations;  these 
configurations  are  correspondingly  labeled  I:  “low  chirp”,  II:  “optimal  chirp”,  III:  “high  chirp”. 

We  then  present  an  analysis  of  the  evolution  of  SNR  across  the  propagation  axis  that  allows  us  to 
identify  changes  in  SNR  which  occur  during  different  stages  of  amplification:  initial  growth , 
saturation ,  and  over-saturation ,  using  a  one-dimensional  OPCPA  model  that  considers  all 
longitudinal  modes  and  associated  noise. 

Task  2  -  Coherent  pulse  synthesis 

Pulse  synthesis  from  two-stage  optical  parametric  amplification 

This  approach  was  followed  in  Milano  and  aimed  at  the  generation  of  a  gap-free  nearly  single¬ 
cycle  pulse  by  coherent  synthesis  of  two  OP  As.  In  the  last  phase  of  the  project,  we  have  developed  the 
first  of  the  schemes  proposed  in  section  “ Methods ,  assumptions  and  Procedures” of  the  third  report, 
based  on  two  separate  OP  As  operating  in  parallel.  Among  the  three  proposed  schemes,  we  selected 
the  approach  based  on  two  parallel  and  independent  OPAs  for  the  amplification  of  an  ultrabroadband 
seed  extending  from  500  to  1000  nm  [Figure  1(a)].  The  OPAs  gain  band  widths  cover  two  partially 
overlapping  spectral  portions  of  this  region,  allowing  the  generation  of  ultrabroadband  pulses  without 
detrimental  gaps  in  the  spectrum  [Figure  1(b)]. 
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Figure  1:  (a)  setup  for  separate  optical  parametric  amplification,  compression  and  coherent 
combination  of  visible  and  near-infrared  pulses.  SHG:  second  harmonic  generation;  WLG:  white 
light  generation,  (b)  Spectra  from  a  typical  OPA I  (visible  NOPA)  and  OPA  II  (degenerate  OPA). 

The  three  building  blocks  of  this  approach  are: 

The  two-stage  IR  OPA,  powered  by  the  fundamental  of  our  laser  source  at  800  nm  and  seeded 
by  the  IR  portion  of  a  WLC  produced  in  sapphire.  The  IR  light  arising  from  this  system  is  then  used  to 
generate  a  supercontinuum  by  self-phase-modulation  in  a  YAG  plate.  The  supercontinuum,  ranging 
from  500  to  1000  nm.;  will  seed  the  two  OP  As  described  in  the  following. 

The  SH-pumped  non-collinear  OPA  (NOPA),  which  amplifies  broadband  visible  pulses  in  the 
500-700  nm  range  [21]  (see  figure  1  (b)); 

The  SH-pumped  OPA  at  degeneracy  [22],  detailed  in  the  First  Interim  Report  of  the  project; 
this  OPA  provides  7-fs  pulses  in  the  700-1000  nm  range  (see  figure  1  (b)). 

The  first  block  aims  at  the  generation  of  the  phase-stable  broadband  white  light;  to  this  end,  a 
great  effort  has  been  devoted  to  the  realization  of  a  reliable  set-up,  driven  by  the  fundamental  beam  of 
our  laser  and  based  on  a  two-stage  IR  OPA.  The  two  OPA  crystals  employ  a  type  II  configuration, 
which  provides  a  stable  narrow-band  amplification.  The  first  OPA  is  tuned  at  1500  nm  and  generates 
the  idler  at  1700  nm.  Since  amplification  arises  from  the  interaction  between  phase-locked  signal  and 
pump  pulses,  the  idler  is  self-phase  stabilized  (23,  24).  The  phase-stable  idler  pulse  is  then  amplified 
in  the  second  stage,  driven  to  saturation  and  designed  in  order  to  optimize  the  idler  energy  stability. 
The  resulting  amplified  beam  is  focused  into  a  YAG  plate  for  the  generation  of  the  ultrabroadband 
supercontinuum.  The  WLC  is  split  in  two  and  used  to  seed  simultaneously  the  visible  NOPA  and  the 
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degenerate  OP  A,  which  are  then  separately  compressed  using  chirped  mirrors  to  nearly  TL  duration. 
The  two  pulses  are  combined  using  a  broadband  beam  splitter  and  their  relative  CEP  and  delay  is 
locked  using  a  nonlinear  correlator. 

Sub-cycle  pulse  synthesis  from  two  OPCPAs 


This  approach  was  followed  at  MIT,  in  collaboration  with  Milano,  and  aimed  at  the  generation  of 
energy  scalable  sub-cycle  optical  waveforms  by  the  coherent  synthesis  of  two  OPCPAs.  Figure  2 
shows  a  schematic  of  the  system.  It  starts  with  an  actively  CEP -stabilized  octave-spanning  Ti:sapphire 
oscillator.  The  oscillator’s  repetition  rate  serves  as  the  master  clock  for  the  full  system,  and  the  1047- 
nm  component  feeds  a  1-kHz  Nd:YLF  chirped  pulse  amplifier  (CPA)  system  to  pump  the  OPCPAs. 
The  oscillator  output  directly  seeds  the  NIR-OPCPA  while  its  spectral  edges  (centered  at  650  nm  and 
930  nm)  undergo  intrapulse  difference-frequency  generation  (DFG)  to  produce  a  pulse  at  2. 15  pm  that 
seeds  the  SWIR-OPCPA.  Using  a  single  oscillator  as  front-end  for  the  entire  system  ensures  the 
coherence  of  the  two  OPCPA  pulses  to  within  environmental  fluctuations  and  drifts  on  subsequent 
beam  paths.  The  designs  of  the  OPCPAs  follow  the  guidelines  described  in  previous  studies  [25,  26] 
for  simultaneously  optimizing  energy  conversion,  amplification  bandwidth,  and  signal-to-noise  ratio. 


Figure  2:  Two  CEP-stabilized,  few-cycle  OPCPAs  centered  at  different  wavelengths  are 
combined  based  on  the  concept  of  coherent  wavelength  multiplexing  to  produce  15-pJ,  0.6-cycle 
pulses  at  1-kHz  repetition  rate.  Full  control  over  the  optical  phase  allows  for  any  optical 
waveform  given  the  amplified  spectrum.  YDFA:  Ytterbium-doped  fiber  amplifier;  BPF: 
bandpass  filter. 


Of  note,  the  inclusion  of  an  acoustooptic  programmable  dispersive  filter  (AOPDF)  in  each  OPCPA 
allows  independent  spectral  phase  and  amplitude  adjustment  of  each  pulse,  enabling  control  and 
optimization  of  the  synthesized  waveform.  Outputs  from  the  two  OPCPAs  are  combined  in  a 
broadband  neutral  beam  splitter.  Besides  the  spectral  phases  (controlled  by  the  AOPDFs),  three  other 
independent  parameters  (see  Fig.  2)  determine  the  synthesized  electric-field  waveform:  the  CEP  of  the 
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NIR-OPCPA  pulse  (cpi),  the  CEP  of  the  SWIR-OPCPA  pulse  (cp2),  and  the  relative  timing  between  the 
two  OPCPA  pulses  (At). 

Precise  stabilization  of  these  three  parameters  is  required  for  coherent  synthesis  of  the  two 
OPCPA  pulses,  and  subsequent  control  of  each  parameter  allows  precise  waveform  shaping.  While 
the  CEP  of  the  SWIR-OPCPA  is  passively  stabilized  due  to  the  intrapulse  DFG  process  used  to 
produce  its  seed,  an  active  feedback  loop  on  the  oscillator  is  implemented  to  ensure  the  CEP  stability 
of  the  NIR-OPCPA.  A  feedback  loop  based  on  a  balanced  cross-correlator  (BCC)  [27]  is  implemented 
to  synchronize  the  two  pulses,  allowing  attosecond-precision  relative  timing  stability.  A  BCC  is  the 
optical  equivalent  of  a  balanced  microwave  phase  detector,  and  is  particularly  suitable  for  timing  drift 
measurements  with  sub-cycle  precision  because  the  balanced  detection  cancels  the  amplitude  noise. 
Once  the  BCC-assisted  feedback  loop  stabilizes  the  relative  timing  between  the  two  OPCPA  pulses,  a 
two-dimensional  spectral-shearing  interferometer  (2DSI)  [28]  is  used  to  measure  the  frequency- 
dependent  group-delay  of  the  synthesized  pulse.  A  2DSI  is  a  variation  of  spectral-shearing 
interferometry,  and  it  circumvents  the  challenge  of  calibrating  interferometer  delay  by  encoding  the 
group-delay  information  in  pure  sinusoidal  fringes  along  a  wavelength  independent  axis,  obtained  by 
scanning  the  relative  phase  of  the  two  spectrally  sheared  components  over  a  few  NIR  periods. 

Results  and  discussion 

Task  1  -  Generation  of  few-optical-cycle  CEP  stable  pulses  from  OPAs 

Simulations  of  noise  evolution  in  OPCPAs 

Integration  of  the  described  nonlinear  equations  allowed  to  follow  the  fluctuating  fields  during 
amplification.  A  graphical  representation  of  the  statistics  of  signal  field  is  given  in  Fig.  3  (a)  and  (b), 
where  we  show  Bm  (vectors),  and  Am  and  nm  (scatter)  for  three  modes  oom  of  the  signal  field.  The 
simulations  refer  to  an  ultra-broadband  OPCPA  system  known  to  be  sensitive  to  PSF  [25,  29,].  The 
amplifier,  pumped  by  a  9-ps  FWHM  Gaussian  pulse  at  1.047  pm  and  seeded  by  a  broadband  (69-THz 
FWHM)  pulse  at  2.094  pm  for  operation  around  degeneracy,  uses  a  3-mm  long  PPSLT  crystal  with 
poling  period  A  =  31.2pm.  These  parameters  are  close  to  the  experimental  conditions  of  Ref.  [25].  In 
this  example,  the  pump-to-seed  energy  ratio  is  106.  The  variances  am2  of  the  noise  field  are 
determined  by  the  quantum  fluctuations  due  to  the  longitudinal  modes  of  the  100-ps-long  simulation 
window.  This  number  is  further  increased  by  a  factor  equal  to  the  number  of  transverse  modes 
amplified  assuming  a  pump  beam  of  100-pm  radius  in  the  3-mm  long  crystal,  which  is  estimated  as 
about  25.  We  note  that  this  choice  results  in  a  calculated  amplified  pulse  contrast  that  closely  matches 
that  observed  in  equivalent  experiments  [25].  For  each  configuration,  we  evaluated  50  independent 
trajectories  triggered  by  uncorrelated  noise  fields.  The  averages  taken  over  this  ensemble  of  classical 
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solutions  correspond  to  quantum-mechanical  expectation  values  [30].  The  results  of  a  batch  of 
simulations  are  depicted  in  Fig.  3(c).  The  incoherent  nature  of  the  amplified  noise  is  evidenced  in 
panel  (d),  representing  the  same  fields  after  compression. 


R 

1° 

R 


Real  part 


Figure  3  (a)  Schematic  representation  of  the  deterministic  part  (arrow)  and  50  stochastic 
components  (circles)  of  a  field  mode,  (b)  Initial  signal  field  distribution  (scatter),  evaluated  at 
three  modes  of  frequency  com,  experiencing  respectively  the  highest  gain  G0,  G0/2  and  G0/10 
and  separated  by  a  phase  shift  imparted  by  the  chirp,  (c)  Depiction  of  the  same  modes  after 
amplification  and  (d)  after  compression;  pedestal  fields  deduced  after  subtracting  the 
deterministic  components  are  indicated  in  the  dashed  circle.  All  data  refer  to  configuration  II. 
Panels  (a)  and  (b)  are  magnified  500  times  with  respect  to  (c)  and  (d). 


The  role  and  amplitude  of  noise  during  amplification  along  propagation  distance  z  was  evaluated 
by  the  two  approaches  introduced  in  previous  reports  and  here  summarized:  in  a  first  approach,  we 
calculated  the  signal-to-noise  ratio  (SNR)  as  the  ratio  between  the  mean  and  the  standard  deviation  of 
the  seed  pulse  energy,  <£'(z)>/A£'(z);  a  second  figure  of  merit  is  the  Signal-to-Pedestal  Ratio  (SPR), 
which  directly  compares  the  average  energy  of  the  incoherent  pedestal  with  the  coherent  pulse.  It  is 
important  to  stress  that  SPR  is  accessible  only  by  numerical  approaches,  since  it  cannot  be 
experimentally  measured  and  the  two  contributions  cannot  be  separated.  Analysis  of  these  figures 
during  amplification  allows  deducing  many  properties  of  the  noise  added  to  the  amplifier  in  the 
presence  of  pump  depletion  and  as  a  function  of  the  amplification  parameters,  such  as  signal  chirp  and 
pump  peak  intensity.  In  particular  we  evaluated  3  configurations,  which  correspond  to:  an  under¬ 
chirped  amplifier,  with  maximum  amplified  signal  bandwidth  but  limited  conversion  efficiency 
(Configuration  I);  an  amplifier  chirped  for  maximum  efficiency-bandwidth  product  (Configuration 
II);  and  an  over-chirped  amplifier  (Configuration  III),  with  excellent  conversion  efficiency  but 
significant  spectral  narrowing.  Mean  and  standard  deviation  of  seed  and  idler  pulse  energies,  together 
with  the  pedestal  energy,  are  given  in  Fig.  4.  The  two  panels  allow  comparing  configurations  I  and  II, 
showing  that  the  pedestal  dramatically  depends  on  amplification  parameters,  which  thus  have  to  be 
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carefully  tuned  in  order  to  minimize  detrimental  tails  in  the  signal  temporal  profile. 


z  [mm] 

Figure  4:  Evolution  of  energy  mean  and  standard  deviation  for  signal  and  idler;  the  energy 
growth  of  the  pedestal  mean  is  also  given.  In  the  case  of  configuration  II  we  indicate  the 
coordinate  at  which  the  pump  peak  is  fully  depleted  (see  inset).  Trends  calculated  for 
configuration  III  (not  shown)  are  comparable  to  the  ones  of  configuration  II 

Evaluation  of  SNR  and  SPR  shows  that  the  two  approaches  for  noise  estimation  follow  different 
behaviors:  in  particular  the  low-chirp  regime  (configuration  I),  exhibits  a  poor  SPR,  due  to  a  strong 
pedestal  component.  This  can  be  explained  by  the  fact  that  energy  fluctuations  A E(z)  and  pedestal 
energy  <EP(z)>  are  not  directly  correlated:  saturation  generally  mitigates  energy  fluctuations,  also 
when  this  energy  mainly  comes  from  noise.  In  this  sense,  the  low-chirp  configuration  is  particularly 
sensitive  because  the  seed  pulse  efficiently  extracts  all  the  energy  from  the  peak  of  the  pump,  but  the 
rest  of  the  pump  energy  is  transferred  to  the  pedestal.  The  effect  of  the  pedestal  is  confirmed  by 
comparing  the  temporal  shape  of  the  pulses  after  compression:  the  uncoherent  pedestal  is  only 
marginally  affected  by  compression,  thus  giving  rise  to  a  strong  plateau  with  the  same  duration  of  the 
pump  beam.  This  long  pulse  has  detrimental  effects  for  those  processes  requiring  a  very  high  pulse 
contrast. 

Task  2  -  Coherent  pulse  synthesis 

Synthesis  from  two-stage  optical  parametric  amplification 

The  first  building  block  of  this  scheme  required  the  preliminary  implementation  of  a  two-stage  IR 
OPA  for  the  generation  of  self-phase  stabilized,  high-energy  pulses;  such  pulses  will  then  drive 
supercontinuum  generation  in  a  thick  YAG  plate.  A  detailed  schematic  of  the  setup  is  shown  in  figure 
5:  a  small  fraction  of  the  driving  pulse  at  800  nm  is  spectrally  broadened  in  a  2-mm  thick  sapphire 
plate;  the  IR  components  of  this  supercontinuum  are  then  amplified  in  an  IR  OPA  based  on  3 -mm 
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thick  BBO  crystal  cut  for  type  II  operation.  We  chose  to  propagate  both  pump  and  seed  with 
extraordinary  (horizontal)  polarization  in  order  to  generate  idler  photons  with  ordinary  (vertical) 
polarization,  as  required  by  the  subsequent  stages.  In  addition,  type  II  was  chosen  because  it  provides 
narrowband  amplification  and  exhibits  higher  gain  thanks  to  the  trapping  effect  induced  by  the 
favourable  group -velocity  mismatches.  In  this  case,  signal  and  idler  pulses  walk  in  opposite  direction 
with  respect  to  the  pump,  so  that  a  nonlinear  interaction  mechanism  localizes  them  under  the  pump 
pulse  and  the  gain  grows  exponentially  even  for  crystal  lengths  well  in  excess  of  the  pulse  splitting 
length.  To  qualitatively  understand  this  trapping  effect,  we  can  consider  the  situation  in  which  the 
signal  pulse  has  moved  slightly  to  the  front  and  the  idler  pulse  to  the  back  of  the  pump  pulse:  during 
the  parametric  process,  the  signal  pulse  generates  idler  photons,  which  move  to  the  back,  i.e.,  toward 
the  peak  of  the  pump;  on  the  other  hand  the  idler  pulse  generates  signal  photons  which  in  turn  move  to 
the  front,  again  toward  the  peak  of  the  pump. 

The  first  OPA  operates  at  1500  nm  signal  wavelength,  corresponding  to  idler  pulses  at  1.7  microns 
wavelength.  Typical  signal+idler  energy  is  of  the  order  of  5  pJ.  After  rejecting  the  residual  pump 
pulse  by  means  of  a  long-pas  filter,  and  the  signal  thanks  to  a  polarizer,  the  idler  pulses  are  sent  to  the 
second  stage  for  further  amplification.  This  stage  employs  a  4-mm  thick  BBO  crystal,  and  is  pumped 
by  800-nm  pulses  with  150  pJ  energy. 


Figure  5:  Detailed  experimental  setup  of  the  two-stage  IR  OPA  for  the  generation  of  gap-free 
supercontinuum  light. 

Signal  and  idler  spectra  from  the  two-stage  OPA  are  shown  in  Fig.  6  (a).  Note  that  here  Idler  and 
Signal  are  named  after  the  first  stage,  and  that  the  second  stage  is  tuned  in  order  to  amplify  the  idler 
light  (black  solid  line).  A  small  non-collinear  angle  between  pump  and  idler  in  the  second  stage 
allows  spatial  separation  of  the  amplifed  idler  from  the  residual  pump  and  signal  photons,  which  may 
be  detrimental  for  the  generation  of  a  phase-stable  supercontinuum.  The  stage  provides  vertically- 
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polarized,  25 -pJ  idler  pulses,  with  energy  fluctuation  of  the  order  of  2%.  Such  pulses  are  then  focused 
in  a  3 -mm  thick  sapphire  plate  for  spectral  broadening  thanks  to  self-phase  modulation.  The  resulting 
supercontinuum  is  shown  as  blue  solid  line  in  Fig.  6(b);  this  spectrum  covers  the  amplification 
bandwidth  of  the  subsequent  visible  NOPA  and  the  degenerate  800-nm  OP  A,  as  shown  by  the  green 
and  red  solid  line  reported  in  the  same  panel. 


Figure  6:  (a)  signal  and  idler  spectra  from  the  second  stage  of  the  double-OPA  setup.  The 
idler  pulses,  with  vertical  polarization,  drive  the  supercontinuum  generation  in  a  4-mm  thick 
YAG  plate,  (b)  Spectrum  of  the  gap-free  supercontinuum  (blue  line)  generated  in  the  YAG 
plate,  compared  to  the  gain  of  the  visible  and  infrared  OPAs  that  will  be  employed  for  its 
amplification. 

The  supercontinum  is  subsequently  split  into  two  identical  replicas  by  a  broadband  ultrathin 
beamsplitter,  and  synchronized  with  the  visible  and  IR  OPAs  for  amplification.  The  visible  NOPA  is 
pumped  by  the  second  harmonic  of  the  800-nm  beam,  provided  by  frequency  doubling  in  a  4-mm 
thick  BBO  crystal.  Such  a  thick  BBO  crystal  was  chosen  in  order  to  narrow  the  SH  spectrum  and 
increase  its  pulse  duration  to  facilitate  its  overlap  with  the  strongly  chirped  seed.  Amplification 
provides  1-pJ  visible  pulses  extending  from  500  nm  to  750  nm  (see  Fig.  7(a),  green  solid  line); 
temporal  compression  is  obtained  by  12  bounces  onto  specially  designed  Double  Chirped  Mirrors 
(DCMs). 

The  degenerate  OPA  has  a  similar  design,  the  only  difference  being  the  angle  between  pump  and 
seed  and  the  required  bounces  on  the  DCMs.  Degenerate  OPA  requires  that  the  seed  propagates 
collinearly  with  the  pump;  here  a  small  non-collinear  angle  allows  to  fulfill  this  condition  and  to 
separate  the  amplified  signal  from  the  spectrally  overlapped  idler.  Compression  is  achieved  by  onto 
DCMs,  designed  to  compensate  dispersion  from  650  nm  to  1.1  microns;  thanks  to  their  capability  to 
introduce  strong  negative  dispersion,  and  to  a  smaller  material  dispersion  in  the  infrared  spectral  range 
than  in  the  visible,  only  two  bounces  are  sufficient  to  obtain  a  duration  close  to  the  TL. 

To  ensure  good  spatial  overlap  of  the  collimated  beams,  the  two  OPAs  are  designed  in  order  to 
employ  the  same  focal  lengths  and  to  propagate  the  beams  for  the  same  distances.  The  two  amplified 
pulses  are  then  synchronized  by  a  delay  line  equipped  with  a  piezoelectric  actuator,  and  collinearly 
combined  by  a  second  ultrathin  beam-splitter.  A  gap-free  spectrum  arising  from  the  combination  of 
the  two  pulses  is  shown  in  Fig.  7(a);  it  supports  sub-4  fs  pulse  duration.  The  pulse  energy  of  the 
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synthesized  pulse  is  500-800  nJ. 


Figure  7:  (a)  Spectrum  of  the  OPAs  light  and  of  the  synthesized  pulses  (blue  line,  log  scale); 
the  total  spectrum  supports  sub-4  fs  pulse  duration,  (b)  Interference  fringes  arising  when  the 
pulses  are  combined  with  200-fs  delay. 


The  last  challenging  step  of  the  pulse  synthesis  is  the  coherent  combination  of  the  two  pulses, 
which  calls  for  careful  control  of  their  relative  delay  and  phase.  Thanks  to  the  spectral  overlap  of  the 
two  beams,  we  could  directly  use  spectral  interferometry  to  characterize  their  delay  fluctuations  (see 
Fig.  7(b));  for  this  purpose,  the  delay  between  the  unlocked  pulses  was  increased  to  about  200  fs  and 
the  spectral  fringes  were  monitored  for  one  minute,  giving  slow  delay  fluctuations  with  rms  of  the 
order  of  3  fs.  When  seeded  with  CEP  stable  continua,  we  plan  to  lock  the  relative  delay  of  the  pulses 
by  using  the  balanced  nonlinear  cross-correlator  (BCC)  [27]  shown  in  Fig.  8,  which  allows 
attosecond-precision  relative  timing  stability  thanks  to  the  capability  of  the  balanced  detection  scheme 
to  cancel  the  amplitude  noise  [31]. 


Figure  8:  Balanced  cross- correlator  (BCC)  used  to  lock  the  delay  of  the  pulses  from  OPA I  and 
OPA  II.  Lij2:  Fused  silica  plates  to  obtain  respectively  compression  and  positive  chirp;  C: 
compensation  plates  to  match  the  dispersion  of  all  arms;  SFG:  sum-frequency  generation 
stages;  PID:  low-pass  filter  and  control  module  for  the  feedback  on  the  delay  line.  The  cartoon 
displays  the  effects  of  a  delay  on  the  OPA  II  pulse. 


14 


Sub-cycle  pulse  synthesis  from  two  OPCPAs 

In  the  following  we  present  the  results  obtained  at  MIT,  in  collaboration  with  Milano,  on  the 
generation  of  sub-cycle  optical  waveforms  by  coherent  synthesis  of  two  OPCPAs.  The  overall 
spectrum  spans  over  1.8  octaves  (green  lines  in  Fig.  9(a))  and  the  energy  of  the  synthesized  pulse  is 


15  pJ.  Due  to  the  gap  in  the  center  of  the  optical  spectrum,  raw  data  of  a  2DSI  measurement  is 


segmented  into  two  parts  and  then  presented.  Figures  9(d)  and  9(e)  demonstrate  the  CEP  stability  of 
the  two  individual  pulses,  with  r.m.s.  fluctuations  as  low  as  135  mrad  and  127  mrad,  respectively. 
Figure  9(f)  characterizes  the  relative  timing  stability.  With  the  feedback  control  of  the  SWIR- 
OPCPA’s  path  length  over  a  bandwidth  of  30  Hz,  the  relative  timing  drift  is  reduced  to  250  as,  less 
than  5%  of  the  oscillation  period  of  the  SWIR-OPCPA  (7.2  fs).  Figures  9(b)  and  9(c)  show  the  raw 
data  of  a  2DSI  measurement  while  Fig.  9(a)  plots  (black  lines)  the  extracted  frequency-dependent 
group-delay  of  the  synthesized  pulse,  which  is  the  derivative  of  the  spectral  phase  with  respect  to 
frequency.  The  2DSI  measurement  shows  that  the  two  OPCPA  pulses  are  temporally  overlapped  and 
each  is  well  compressed  to  within  10%  of  its  transform-limited  pulse  duration.  In  our  system,  the 
CEPs  can  be  varied  by  slight  tuning  of  any  dispersive  element,  including  the  AOPDFs.  The  values  of 
the  CEP  will  be  determined  automatically  in  situ  when  strong-field  experiments  are  conducted  and 
hence  CEP  tunability  is  sufficient  from  an  experimental  point  of  view. 
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Figure  9  (a)  Optical  spectrum  and  frequency-dependent  group-delay  of  the  synthesized  pulses. 
The  overall  spectrum  spans  over  1.8  octaves  and  supports  sub-cycle  pulses,  (b)  is  the  2DSI 
trace  for  the  NIR-OPCPA  and  (c)  is  that  for  the  SWIR-OPCPA.  CEP  stabilities  are  verified 
using  nonlinear  interferograms.  (d)  f-2f  interferogram,  measuring  135  mrad  rms  (5 -shot 
integration)  CEP  fluctuations  for  the  NIR-OPCPA.  (e)  f-3f  interferogram,  measuring  127  mrad 
rms  (5 -shot  integration)  CEP  fluctuations  for  the  SWIR-OPCPA.  The  BCC-assisted  feedback 
loop  guarantees  the  relative  timing  stability  and  (f)  shows  BCC  measurements  of  the  free- 
running  system  (black)  and  the  closed-loop  system  (red).  The  closed-loop  system  ensures  a 
relative  timing  drift  of  250  as.  The  relative  timing  drift  could  be  reduced  even  further  to  100  as 
if  the  feedback  bandwidth  were  extended  to  100  Hz. 
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Figure  10  plots  a  synthesized  electric-field  waveform  and  intensity  profile  assuming  the  CEPs 
(cpi=650mrad,  (p2=-750mrad)  optimal  for  achieving  the  shortest  high  field  transient,  which  lasts  only 
0.8  cycles  (amplitude  FWHM)  of  the  carrier  (centroid)  frequency  (kc  =  1.26  pm).  The  lower  inset  of 
Fig.  10(a)  clearly  shows  that  the  synthesized  electric-field  waveform  is  non-sinusoidal  and  the  main 
feature  lasts  less  than  an  optical  cycle.  As  an  example  of  waveform  shaping  made  possible  by  tuning 
parameters  of  our  system,  Figs.  10(b)  and  10(c)  show  two  waveforms  as  the  CEP  and  the  relative 
timing  are  changed.  Due  to  the  large  gap  in  the  combined  spectrum,  there  are  wings  4.8  fs  from  the 
central  peak  as  shown  in  Fig.  10(a).  These  wings  should  be  absent  in  the  approach  follwed  in  Milano, 
since  the  synthesized  pulses  are  spectrally  overlapped,  thus  generating  a  gap-free  spectrum.  However, 
for  processes  initiated  by  strong-field  ionization,  these  wings  have  a  negligible  effect.  For  more 
demanding  applications,  the  wings  can  be  suppressed  by  extension  of  the  coherent  wavelength 
multiplexing  scheme  to  include  a  third  OPCPA,  centered  at  1.5  pm,  to  fill  the  spectral  gap.  The 
synthesized  waveforms  are  important  for  optimizing  the  HHG  process,  which  is  to  date  the  only 
demonstrated  technique  for  generating  isolated  attosecond  pulses  [9]. 


Figure  10:  (a)  Electric-field  waveform  of  the  synthesized  sub-cycle  pulses,  calculated 
assuming  CEPs  and  a  relative  timing  of  cpi=650  mrad,  cp2=-750  mrad,  At=0.0  fs  optimal  for 
achieving  the  shortest  high-field  transient.  Lower  inset:  the  waveform  is  superimposed  with  the 
electric  field  oscillating  at  the  carrier  frequency:  the  synthesized  electric-field  waveform  lasts 
less  than  an  optical  cycle.  Upper  inset:  corresponding  intensity  profile,  (b,  c)  Waveforms  under 
r.m.s.  residual  jitters.  While  the  red  solid  line  is  the  unperturbed  waveform  of  panel  (a),  the 
black  dotted  line  is  obtained  by  adding  127  mrad  to  cp2  in  (b),  and  250  as  to  At  in  (c). 

As  an  example,  we  numerically  solve  the  time-dependent  Schrodinger  equation  (TDSE)  for  a  He- 
atom  in  a  strong  laser  field  to  illustrate  a  possible  use  of  our  source  for  driving  direct  isolated  soft-x- 
ray  pulse  generation.  The  achievable  peak  intensity  (6xl014  W/cm2)  is  chosen  such  that  the  total 
ionization  is  below  the  critical  ionization  level  in  helium.  With  choice  of  CEPs  and  timing  as  in  Fig. 
10(a),  substantial  ionization  is  limited  to  one  optical  half-cycle  and  an  isolated  soft-x-ray  pulse 
spanning  over  250  eV  is  generated  without  the  need  for  gating  techniques  [11]  or  spectral  filtering 
which  typically  limit  the  obtainable  bandwidth.  Using  an  additional  Sn  filter,  which  blocks  the  strong 
IR  driving  field  and  the  nonlinearly  chirped  low-photon-energy  spectral  content  below  70  eV,  leads  to 
an  isolated  150-as  pulse  centered  at  200  eV.  Of  note,  the  non-sinusoidal  electric-field  waveform  leads 
to  drastically  changed  electron  trajectories  (compared  to  those  from  a  sinusoidal  driving  field) 
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resulting  in  corresponding  changes  in  quantum  diffusion  and  atto-chirp,  which  can  be  controlled  by 
means  of  the  waveform  synthesis  parameters  (cpi,  cp2,  and  At).  In  this  example,  quantum  diffusion 
dominates  over  ionization  rate  and  effectively  eliminates  the  radiation  from  long  trajectories,  resulting 
in  isolated  soft-x-ray  pulse  generation  solely  from  short  trajectories.  This  gives  an  example  of  the 
capability  of  our  sub-cycle  waveform  to  simultaneously  isolate  the  ionization  process  and  manipulate 
electron  trajectories  within  an  optical  cycle,  allowing  unprecedented  control  of  the  HHG  process. 
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Conclusions 


In  the  fourth  and  final  reporting  period  of  this  project  we  have  concluded  the  theoretical  studies 
on  the  effects  of  quantum  noise  on  OPCPAs,  and  optimized  two  experimental  set-ups  for  coherent 
pulse  synthesis  and  CEP-stable  single-cycle  pulse  generation,  developed  in  Milano  and  at  MIT 
respectively. 

In  Task  1  we  have  analyzed  noise  evolution  in  OPCPAs,  obtaining  useful  guidelines  for  their 
design  which  is  a  preliminary  step  to  the  pulse  synthesis.  We  have  performed  a  quantum-mechanically 
consistent  numerical  investigation  of  the  dynamics  of  superfluorescence  growth  in  a  realistic  high- 
gain  OPCPA.  Thanks  to  the  model’s  capability  to  simulate  also  the  saturation  stage  of  amplification, 
this  investigation  for  the  first  time  captured  all  dynamics  of  a  quantum-noise-contaminated  OPCPA, 
addressing  the  important  practical  issues  of  signal  energy  stability  and  pulse  contrast.  These  quantities 
display  different  evolution  dynamics  throughout  the  amplification  process.  Three  operating  conditions 
were  explored,  characterized  by  different  chirps  of  the  input  seed.  We  find  that  the  chirp  maximizing 
the  efficiency-bandwidth  product  is  also  characterized  by  the  smallest  contribution  of  the  noise. 
Significantly,  we  find  that  while  amplifier  saturation  improves  the  signal’s  shot-to-shot  energy 
stability,  it  does  not  necessarily  improve  the  pulse  contrast.  Knowledge  of  these  dynamics  increases 
our  fundamental  understanding  of  quantum  noise  in  parametric  amplification;  it  also  provides 
important  insight  for  the  optimization  of  OPCPA  systems  applied  to  the  study  of  strong-field  laser 
physics.  These  results  have  been  published  in  the  following  paper: 

Cristian  Manzoni,  Jeffrey  Moses,  Franz  X.  Kartner,  and  Giulio  Cerullo,  ((Excess  quantum  noise 
in  optical  parametric  chirped-pulse  amplification,  ”  Opt.  Express  19,  8357-8366  (2011) 

The  full  text  is  provided  in  the  Appendix. 

In  addition  a  paper  has  been  presented  at  the  Conference  on  Lasers  and  Electro-Optics  (CLEO) 
(Munich,  Germany,  May  22-26  2011)  with  the  title 

J.  Moses,  C.  Manzoni,  G.  Cerullo,  and  F.X.  Kartner,  Superfluorescence  Dynamics  of  OPCPAs 
in  the  Saturation  Regime  ” 

The  abstract  is  provided  in  the  Appendix. 

In  Task  2  we  have  implemented  two  different  experimental  schemes  for  coherent  pulse  synthesis. 
The  results  can  be  summarized  as  follows: 

Synthesis  from  two-stage  optical  parametric  amplification 

We  have  realized  a  two-stage  IR-OPA  for  the  generation  of  a  gap-free  phase-stable 
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ultrabroadband  WLC.  This  setup  removes  the  spectral  gap  at  800  nm  typically  arising  when  spectral 
broadening  the  fundamental  frequency  of  a  Ti:sapphire  laser.  The  WLC  is  split  in  two  and  used  to 
seed  simultaneously  the  visible  NOPA  and  the  degenerate  OPA,  which  are  then  separately 
compressed  using  chirped  mirrors  to  nearly  TL  duration.  The  two  pulses  are  combined  using  a 
broadband  beam  splitter  and  we  plan  to  lock  their  relative  CEP  and  delay  using  a  nonlinear  cross¬ 
correlator.  Preliminary  results  demonstrated  that  combination  of  the  OPA  pulses  allows  obtaining  a 
broadband  pulse  extending  from  550  to  1000  nm,  corresponding  to  a  transform- limited  duration  of  4 
fs.  Spectral  interferometry  demonstrated  that  the  delay  of  the  two  pulses  exhibits  fluctuations  which 
can  be  compensated  by  a  nonlinear  correlator.  The  results  of  these  preliminary  experiments  have  been 
submitted  as  a  contribution  to  the  “Ultrafast  Optics  2011”  conference  (September  26-30,  2011, 
Monterey,  CA)  with  the  title: 

C.  Manzoni,  S.W  Huang,  G.  Cirmi,  J.  Moses,  F.  X.  Kartner,  and  G.  Cerullo,  “Ultrabroadband 
pulse  generation  by  coherent  synthesis  of  two  optical  parametric  amplifiers  ” 

The  abstract  of  the  contribution  is  given  in  the  appendix. 

Sub-cycle  pulse  synthesis  from  two  OPCPAs 

We  have  developed  a  scalable  optical  waveform  synthesyzer  scheme  based  on  fully  controlled 
coherent  wavelength  multiplexing  of  high-energy,  few-cycle  optical  pulses  from  multi-color 
OPCPAs.  Currently,  the  system  generates  a  sub-cycle  waveform  with  a  spectrum  spanning  close  to 
two  octaves  and  15-pJ  pulse  energy.  It  can  be  readily  scaled  both  in  energy  and  bandwidth  given  the 
proven  wavelength  tunability  of  OPCPAs.  A  numerical  study  shows  the  uniqueness  of  our  source  for 
direct  isolated  soft-x-ray  pulse  generation  based  on  HHG,  eliminating  the  need  for  gating  techniques 
or  spectral  filtering.  The  system  is  capable  of  stabilizing  and  controlling  all  independent  parameters 
that  define  the  synthesized  electric-field  waveform;  this  new  high-intensity  laser  architecture  can  be 
applied  to  optical  field-emission,  tunneling  ionization  studies,  time -resolved  spectroscopy,  and  in 
general,  attosecond  control  of  strong-field  physics  experiments.  The  results  of  this  work  have  been 
published  in  the  following  paper: 

S.-W.  Huang,  G.  Cirmi,  J.  Moses,  K.-H.  Hong,  S.  Bhardwaj,  J.  R.  Birge,  L.-J.  Chen,  E.  Li,  B.  J. 
Eggleton,  G.  Cerullo,  and  F.  X  Kartner,  “Scalable  High-Energy  Sub-Cycle  Wayeform  Synthesis  for 
Strong-Field  Physics  ”,  accepted  for  publication  in  Nature  Photonics  and  currently  available  on-line. 
The  full  text  is  provided  in  the  Appendix. 
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High-energy  pulse  synthesis  with  sub-cycle 
waveform  control  for  strong-field  physics 

Shu-Wei  Huang1,  Giovanni  Cirmi1,  Jeffrey  Moses1,  Kyung-Han  Hong1,  Siddharth  Bhardwaj1, 
Jonathan  R.  Birge1,  Li-Jin  Chen1,  Enbang  Li2,  Benjamin  J.  Eggleton2,  Giulio  Cerullo3 
and  Franz  X.  Kartner1-4 ** 


Over  the  last  decade,  control  of  atomic-scale  electronic  motion 
by  non-perturbative  optical  fields  has  broken  tremendous  new 
ground  with  the  advent  of  phase-controlled  high-energy  few- 
cycle  pulse  sources1.  The  development  of  close  to  single¬ 
cycle,  carrier-envelope  phase  controlled,  high-energy  optical 
pulses  has  already  led  to  isolated  attosecond  EUV  pulse  gener¬ 
ation2,  expanding  ultrafast  spectroscopy  to  attosecond  resol¬ 
ution1.  However,  further  investigation  and  control  of  these 
physical  processes  requires  sub-cycle  waveform  shaping, 
which  has  not  been  achievable  to  date.  Here,  we  present  a 
light  source,  using  coherent  wavelength  multiplexing,  that 
enables  sub-cycle  waveform  shaping  with  a  two-octave-spa n- 
ning  spectrum  and  a  pulse  energy  of  15  |mJ.  It  offers  full 
phase  control  and  allows  generation  of  any  optical  waveform 
supported  by  the  amplified  spectrum.  Both  energy  and  band¬ 
width  scale  linearly  with  the  number  of  sub-modules,  so  the 
peak  power  scales  quadratically.  The  demonstrated  system  is 
the  prototype  of  a  class  of  novel  optical  tools  for  attosecond 
control  of  strong-field  physics  experiments. 

Since  the  invention  of  pulsed  lasers,  the  ultrafast  laser  science 
community  has  strived  for  ever  broader  optical  bandwidths, 
shorter  pulse  durations,  higher  pulse  energies  and  improved  phase 
control.  Each  breakthrough  in  generation  methods  has  led  to  new 
scientific  discoveries  in  a  wide  range  of  fields3-5.  Recent  investi¬ 
gations  of  phenomena  at  the  intersection  of  ultrafast  and  strong- 
field  laser  physics,  such  as  high-harmonic  generation  (HHG)6  and 
strong- field  ionization7,  have  demanded  that  laser  sources 
combine  each  of  the  breakthroughs  mentioned  above. 
Investigation  and  control  of  the  strong- field  light-matter  interaction 
simultaneously  requires  a  multi-octave-spanning  bandwidth,  an  iso¬ 
lated  sub-cycle  waveform,  peak  intensity  above  1  x  10 14  W  cm-2 
and  full  phase  control.  Such  features  would  allow  arbitrary 
shaping  of  the  strong  electric-field  waveform  for  steering  ionized 
electron  wave  packets8  and  precise  control  of  tunnelling  and  multi¬ 
photon  ionization  events. 

For  over  two  decades,  laser  scientists  have  sought  to  extend  laser 
bandwidths  and  achieve  sub -cycle  optical  waveforms  by  synthesiz¬ 
ing  multiple  laser  sources9.  Attempts  to  combine  two  independent 
mode-locked  lasers  have  met  with  some  success,  for  example  in  fre¬ 
quency  metrology10,11,  but  are  challenging  because  of  the  differential 
phase  noise  beyond  the  achievable  feedback  loop  bandwidth.  This 
problem  was  recently  circumvented  by  coherently  adding  two 
pulse  trains  derived  from  the  same  fibre  laser,  resulting  in  the  first 
demonstration  of  an  isolated  single-cycle  optical  pulse  source12. 
This  proved  the  feasibility  of  pulse  synthesis  at  the  nanojoule 


level,  but  achieving  high  pulse  energy  requires  the  synthesis  of 
low-repetition-rate  pulses,  which  is  a  challenge  because  of  the 
environmental  perturbations  typical  of  high-energy  amplifiers.  An 
approach  to  high-energy  pulse  synthesis  based  on  combining  the 
pump,  signal  and  idler  of  a  multi- cycle  optical  parametric  amplifier 
is  being  investigated,  and  shows  the  potential  to  produce  multiple 
single-cycle  pulses  under  a  multi-cycle  envelope  with  pulse  separ¬ 
ation  on  the  order  of  a  few  femtoseconds13. 

In  this  Letter,  we  address  the  challenge  of  high-energy  sub-cycle 
optical  waveform  synthesis.  We  demonstrate  a  new  approach,  based 
on  coherent  wavelength  multiplexing  of  ultra-broadband  optical 
parametric  chirped  pulse  amplifiers  (OPCPAs),  for  the  generation 
of  fully  controlled  high-energy  non-sinusoidal  optical  waveforms 
with  spectra  spanning  close  to  two  octaves.  By  means  of  simulation, 
we  present  an  example  of  the  unique  features  of  our  source  as  a 
driver  for  isolated  strong-field  physics  experiments:  the  confinement 
of  the  strong-field  light-matter  interaction  to  within  an  optical  cycle 
and  attosecond  control  of  the  interaction.  The  system  coherently 
combines  two  carrier-envelope  phase  (CEP) -controlled,  few-cycle 
pulses  obtained  from  different  OPCPAs:  (i)  a  near-infrared  (NIR) 
OPCPA,  producing  25  pj,  9  fs  pulses  centred  at  870  nm  and  (ii)  a 
short-wavelength  infrared  (SWIR)  OPCPA,  producing  25  pj,  24  fs 
pulses  centred  at  2.15  pm.  The  ultra-broadband  OPCPA  is  the 
most  promising  technology  for  producing  wavelength-tunable, 
high-peak-power  and  high-average-power,  few-cycle  optical  pulses 
with  good  pre-pulse  contrast14.  Furthermore,  an  ultra-broadband 
OPCPA  maintains  good  CEP  stability  due  to  the  low  thermal  load 
and  the  small  dispersion  required  to  stretch  and  compress 
the  signals. 

Figure  1  shows  a  schematic  of  the  system.  It  starts  with  an  actively 
CEP -stabilized  octave-spanning  Ti:sapphire  oscillator.  Using  a 
single  oscillator  as  the  front-end  for  the  entire  system  ensures  the 
coherence  of  the  two  OPCPA  pulses  to  within  environmental  fluc¬ 
tuations  and  drifts  on  subsequent  beam  paths.  The  designs  of  the 
OPCPAs  follow  the  guidelines  described  in  previous  studies15,16 
for  simultaneously  optimizing  energy  conversion,  amplification 
bandwidth  and  signal-to-noise  ratio.  Of  note,  the  inclusion  of  an 
acousto-optic  programmable  dispersive  filter  (AOPDF)  in  each 
OPCPA  allows  independent  spectral  phase  and  amplitude  adjust¬ 
ment  of  each  pulse,  enabling  control  and  optimization  of  the 
synthesized  waveform. 

Outputs  from  the  two  OPCPAs  are  combined  in  a  broadband 
neutral  beamsplitter.  The  overall  spectrum  spans  over  1.8  octaves 
(green  lines  in  Fig.  2a),  and  the  energy  of  the  synthesized  pulse  is 
15  pj.  Besides  the  spectral  phases  (controlled  by  the  AOPDFs), 
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Figure  1  |  Schematic  of  the  high-energy  optical  waveform  synthesizer.  Two  CEP-stabilized,  few-cycle  OPCPAs  centred  at  different  wavelengths  are  combined 
based  on  the  concept  of  coherent  wavelength  multiplexing  to  produce  a  fully  controlled  non-sinusoidal  optical  waveform  with  a  pulse  energy  of  15  juuJ  at  a 
repetition  rate  of  1  kHz.  Full  control  over  the  optical  phase  allows  for  any  optical  waveform  given  the  amplified  spectrum.  YDFA,  ytterbium-doped  fibre 
amplifier;  BPF,  bandpass  filter;  DFG,  difference  frequency  generation  (intra-pulse). 
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Figure  2  |  Characterization  of  the  synthesized  pulses,  a,  Optical  spectrum  (green)  and  frequency-dependent  group  delay  (black)  of  the  synthesized  pulses. 
The  overall  spectrum  spans  over  1.8  octaves  and  supports  non-sinusoidal  waveforms  with  sub-cycle  features.  b,c,  2DSI  trace  for  the  NIR  OPCPA  (b)  and  the 
SWIR  OPCPA  (c).  The  2DSI  measurements  show  that  the  two  pulses  are  temporally  overlapped  and  well  compressed  to  within  10%  of  the  transform-limited 
pulse  duration.  CEP  stabilities  are  verified  using  nonlinear  interferograms  with  five-shot  integration,  d,  f-2f  interferogram,  measuring  135  mrad  r.m.s.  CEP 
fluctuations  over  30  s  for  the  NIR  OPCPA.  e,  f-3f  interferogram,  measuring  127  mrad  r.m.s.  CEP  fluctuations  over  30  s  for  the  SWIR  OPCPA.  Spatial  properties 
are  characterized  by  measuring  the  beam  profiles  and  the  A/I2  values.  f,g,  Beam  profile  of  the  NIR  OPCPA  (f)  and  the  SWIR  OPCPA  (g).  The  A/I2  value  of  the 
NIR-OPCPA  is  1.2  and  that  of  the  SWIR  OPCPA  is  1.3.  The  BCC-assisted  feedback  loop  guarantees  the  relative  timing  stability,  h,  BCC  measurements  of  the 
free-running  (black)  and  closed-loop  (red)  systems.  The  closed-loop  system  ensures  a  relative  timing  drift  of  250  as,  less  than  5%  of  the  oscillation  period  of 
the  SWIR  OPCPA  (over  10  s). 
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Figure  3  |  The  synthesized  electric-field  waveforms,  a,  Here,  we  assume  CEPs  (<£,  =  650  mrad,  </>2  =  -750  mrad)  optimal  for  achieving  the  shortest 
high-field  transient,  which  lasts  only  0.8  cycles  (amplitude  FWHM)  of  the  carrier  (centroid)  frequency.  Lower  inset:  the  waveform  is  plotted  in  a  shorter  time 
window  and  superimposed  with  an  electric  field  oscillating  at  the  carrier  (centroid)  frequency,  showing  that  the  synthesized  electric-field  waveform  is  non- 
sinusoidal  and  the  main  feature  lasts  less  than  an  optical  cycle.  Upper  inset:  corresponding  intensity  profile.  About  one-third  of  the  pulse  energy  is  contained 
in  the  main  pulse.  As  an  example  of  waveform  shaping  made  possible  by  tuning  the  parameters  of  our  system,  two  additional  atypical  waveforms  are  shown, 
b,  A  waveform  synthesized  by  adding  500  mrad  to  both  ^  and  </>2.  c,  A  waveform  synthesized  by  adding  1  fs  to  At. 


three  other  independent  parameters  (Fig.  1)  determine  the  syn¬ 
thesized  electric-field  waveform:  the  CEP  of  the  NIR  OPCPA 
pulse  (4>i),  the  CEP  of  the  SWIR  OPCPA  pulse  (02)  and  the  relative 
timing  between  the  two  OPCPA  pulses  (A t).  Precise  stabilization  of 
these  three  parameters  is  required  for  coherent  synthesis  of  the  two 
OPCPA  pulses,  and  subsequent  control  of  each  parameter  allows 
precise  waveform  shaping.  Although  the  CEP  of  the  SWIR 
OPCPA  is  passively  stabilized  due  to  the  intrapulse  difference- 
frequency  generation  (DFG)17  used  to  produce  its  seed,  an  active 
feedback  loop  on  the  oscillator  is  implemented  to  ensure  the  CEP 
stability  of  the  NIR  OPCPA.  Figure  2d,e  demonstrates  the  CEP  stab¬ 
ility  of  the  two  individual  pulses,  with  r.m.s.  fluctuations  as  low  as 
135  mrad  and  127  mrad,  respectively.  Figure  2h  characterizes  the 
relative  timing  stability.  A  feedback  loop  based  on  a  balanced  cross¬ 
correlator  (BCC)18  is  implemented  to  synchronize  the  two  pulses, 
allowing  attosecond-precision  relative  timing  stability.  With  the  feed¬ 
back  control  of  the  SWIR  OPCPA’s  path  length  over  a  bandwidth  of 
30  Hz,  the  relative  timing  drift  is  reduced  to  250  as,  less  than  5%  of 
the  oscillation  period  of  the  SWIR  OPCPA  (7.2  fs). 

Once  the  BCC- assisted  feedback  loop  stabilizes  the  relative 
timing  between  the  two  OPCPA  pulses,  a  two-dimensional 
spectral- shearing  interferometer  (2DSI)19  is  used  to  measure  the 
frequency- dependent  group  delay  of  the  synthesized  pulse. 
Figure  2b, c  presents  the  raw  data  of  a  2DSI  measurement,  and 
Fig.  2a  plots  (black  lines)  the  extracted  frequency- dependent 
group  delay  of  the  synthesized  pulse,  which  is  the  derivative  of  the 
spectral  phase  with  respect  to  frequency.  The  2DSI  measurement 
shows  that  the  two  OPCPA  pulses  are  temporally  overlapped,  and 
each  is  well  compressed  to  within  10%  of  its  transform-limited 
pulse  duration. 

In  our  system,  the  CEPs  can  be  varied  by  slight  tuning  of  any  dis¬ 
persive  element,  including  the  AOPDFs20.  The  values  of  the  CEP 
will  be  determined  automatically  in  situ  when  strong-field  exper¬ 
iments  are  conducted21,  so  CEP  tunability  is  sufficient  from  an 
experimental  point  of  view.  In  summary,  our  system  is  capable  of 
stabilizing  and  controlling  all  independent  parameters  that  define 
the  synthesized  electric-field  waveform.  Figure  3a  plots  a  syn¬ 
thesized  electric-field  waveform  and  intensity  profile  assuming  the 
CEPs  (^  =  650  mrad,  (j>2= —75 0  mrad)  optimal  for  achieving 


the  shortest  high-field  transient,  which  lasts  only  0.8  cycles  (ampli¬ 
tude  FWHM)  of  the  carrier  (centroid)  frequency  (Ac=  1.26  pan). 
The  lower  inset  of  Fig.  3a  clearly  shows  that  the  synthesized  elec¬ 
tric-field  waveform  is  non- sinusoidal,  and  the  main  feature  lasts 
less  than  an  optical  cycle.  As  an  example  of  waveform  shaping 
made  possible  by  tuning  the  parameters  of  our  system,  Fig.  3b, c 
shows  two  atypical  waveforms  as  the  CEP  and  relative  timing  are 
changed.  Because  of  the  large  gap  in  the  combined  spectrum, 
there  are  wings  4.8  fs  from  the  central  peak,  as  shown  in  Fig.  3a. 
As  we  will  show  below,  for  processes  initiated  by  strong-field  ioniz¬ 
ation,  these  wings  have  a  negligible  effect.  For  more  demanding 
applications,  the  wings  can  be  suppressed  by  extension  of  the  coher¬ 
ent  wavelength  multiplexing  scheme  to  include  a  third  OPCPA, 
centred  at  1.5  pan  (ref.  22),  to  fill  the  spectral  gap.  The  synthesized 
waveforms  are  important  for  optimizing  the  HHG  process6,  which 
is,  to  date,  the  only  demonstrated  technique  for  generating  isolated 
attosecond  pulses2.  As  an  example,  we  numerically  solve  the  time- 
dependent  Schrodinger  equation  (TDSE)  for  a  helium  atom  in  a 
strong  laser  field  to  illustrate  a  possible  use  of  our  source  for 
driving  direct  isolated  soft  X-ray  pulse  generation  (Fig.  4).  The 
achievable  peak  intensity  (6  x  10 14  W  cm-2)  is  chosen  such  that 
the  total  ionization  is  below  the  critical  ionization  level  in 
helium23.  With  the  choice  of  CEPs  as  in  Fig.  4a,  substantial  ioniz¬ 
ation  is  limited  to  one  optical  half- cycle,  and  an  isolated  soft 
X-ray  pulse  spanning  over  250  eV  is  generated  (Fig.  4b, c)  without 
the  need  for  gating  techniques24  or  spectral  filtering,  which  typically 
limit  the  obtainable  bandwidth.  Using  an  additional  tin  filter,  which 
blocks  the  strong  IR  driving  field  and  the  nonlinearly  chirped  low- 
photon- energy  spectral  content  below  70  eV,  leads  to  an  isolated 
150  as  pulse  centred  at  200  eV.  Of  note,  the  non-sinusoidal  elec¬ 
tric-field  waveform  leads  to  drastically  changed  electron  trajectories 
(compared  to  those  from  a  sinusoidal  driving  field),  resulting  in 
corresponding  changes  in  quantum  diffusion  and  atto- chirp, 
which  can  be  controlled  by  means  of  the  waveform  synthesis  par¬ 
ameters  ((/>!,  4>2  and  A t).  In  this  example,  quantum  diffusion  dom¬ 
inates  over  ionization  rate  (see  Supplementary  Information)  and 
effectively  eliminates  the  radiation  from  long  trajectories,  resulting 
in  isolated  soft  X-ray  pulse  generation  solely  from  short  trajectories 
(Fig.  4b).  This  gives  an  example  of  the  capability  of  our  light  source 
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Figure  4  |  Extreme  nonlinear  optics  with  sub-cycle  manipulated  waveforms.  TDSE  simulation  results  of  the  single-atom  HHG  show  the  uniqueness  of  our 
source  for  direct  isolated  soft  X-ray  pulse  generation,  a,  Ionization  dynamics  (red)  induced  in  helium  by  a  linearly  polarized  electric-field  waveform  (black)  assuming 
a  peak  intensity  of  6  x  1014  W  cm-2,  </>-,  =  960  mrad  and  4>2  =  -440  mrad.  b,  Spectrogram  of  the  HHG  superimposed  with  the  calculated  classical  trajectories. 
Returning  trajectories  from  three  ionization  events  (2,  main  pulse;  1  and  3,  satellite  pulses)  are  shown  for  clear  interpretation  of  the  spectrogram.  The  synthesized 
pulse  isolates  the  ionization  process  to  a  half  optical  cycle,  and  a  continuum  spectrum  spanning  more  than  250  eV  can  be  achieved.  The  isolated  soft  X-ray  pulse 
has  the  same  sign  of  chirp  over  80%  of  the  spectrum,  so  the  compression  setup  can  be  simplified,  c,  Isolated  soft  X-ray  pulse  plotted  in  the  time  domain  before 
(pink)  and  after  (black  line)  a  100-nm-thick  Sn  filter.  The  Sn  filter  is  chosen  for  its  ability  to  block  the  strong  IR  driving  field  and  the  nonlinearly  chirped  low-photon- 
energy  spectral  content,  and  its  good  transmission  in  the  soft  X-ray  range.  The  filtered  isolated  soft  X-ray  pulse  has  a  FWHM  duration  of  150  as. 


to  simultaneously  isolate  the  ionization  process  and  manipulate 
electron  trajectories  within  an  optical  cycle,  allowing  unprecedented 
control  of  the  HHG  process. 

In  conclusion,  we  have  presented  a  scalable  waveform  synthesis 
scheme  based  on  fully  controlled  coherent  wavelength  multiplexing 
of  high-energy,  few-cycle  optical  pulses  from  multi-colour  OPCPAs. 
Currently,  the  system  generates  a  non-sinusoidal  waveform  that  can 
be  used  to  drive  isolated  strong-field  physics  experiments.  The  pulse 
energy  is  15  pj,  with  the  spectrum  spanning  close  to  two  octaves, 
and  it  can  be  readily  scaled  both  in  energy  and  bandwidth  given 
the  proven  wavelength  tunability  of  OPCPAs25  (see 
Supplementary  Information).  A  numerical  study  shows  the  unique¬ 
ness  of  our  source  for  direct  isolated  soft  X-ray  pulse  generation 
based  on  HHG,  eliminating  the  need  for  gating  techniques24  or 
spectral  filtering.  In  addition  to  this  application,  this  new  high- 
intensity  laser  architecture  can  be  applied  to  optical  field- emission26, 
tunnelling  ionization  studies27,  time-resolved  spectroscopy28  and,  in 
general,  attosecond  control  of  strong-field  physics  experiments. 

Methods 

OPCPA  setup.  The  system  schematic  is  presented  in  Supplementary  Fig.  SI. 

Both  OPCPAs  are  pumped  by  an  optically  synchronized  (injection  seeded  by  the 
Octavius-85M  Ti:sapphire  oscillator  from  IdestaQE,  Inc.)  Nd:YLF  chirped  pulse 
amplifier  (CPA),  which  generates  3.5  mj,  12  ps  pulses  at  1,047  nm.  The  SWIR 
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OPCPA,  pumped  by  1  mj  of  the  Nd:YLF  CPA  output,  follows  the  design  method  in 
ref.  15.  The  seed,  produced  by  intrapulse  DFG  of  the  oscillator,  is  first  stretched  by  a 
bulk  silicon  block  to  5  ps  and  pre-amplified  to  1.5  pj  in  the  first  OPCPA  stage 
using  periodically  poled  lithium  niobate  (PPLN).  The  pre-amplified  pulse  is  further 
stretched  to  9.5  ps  by  an  infrared  AOPDF,  amplified  to  25  pj  in  periodically 
poled  stoichiometric  lithium  tantalate  (PPSLT),  and  then  compressed  to  24  fs  in  a 
broadband  anti- reflection  coated  quartz  glass  block  (Suprasil  300).  For  the  NIR 
OPCPA,  a  2  mj  fraction  of  the  Nd:YLF  CPA  output  is  frequency- doubled  in  a 
lithium  triborate  (LBO)  crystal  and  used  to  amplify  the  oscillator  output.  The  seed  is 
first  stretched  to  5  ps  by  a  Brewster  prism  stretcher.  The  signal  is  pre-amplified  in  a 
double-pass  configuration  in  a  type-I,  5-mm-long  [3-barium  borate  (BBO)  crystal  to 
2  pj.  The  amplified  pulse  is  further  stretched  to  6.2  ps  by  an  AOPDF  and  a  grating 
stretcher,  amplified  to  25  pj  in  BBO  and  then  compressed  to  9  fs  in  a  Brewster-cut 
N-LaSF9  block. 

Beam  combining.  The  outputs  of  the  two  OPCPAs  are  combined  in  a  broadband 
neutral  beamsplitter,  which  introduces  25%  energy  loss.  In  addition,  only  half  of 
the  synthesized  pulse  energy  is  available  for  experiments,  because  the  other  half  is 
directed  to  the  balanced  cross-correlator  (BCC).  Because  the  pulse  energy  directed 
to  the  BCC  is  much  greater  than  is  needed,  a  custom-made  dichroic  mirror  can  be 
implemented  to  improve  the  experimentally  available  pulse  energy  by  a  factor  of  2.5. 
Thus,  in  an  optimized  system,  waveform  synthesis  could  be  achieved  with  very 
low  losses. 

We  chose  to  combine  the  two  OPCPA  pulses  in  a  ‘constant  waist  width’ 
fashion29,  which  is  inherently  compatible  with  OPCPA  configurations.  As  shown 
theoretically  in  ref.  29,  the  ‘constant  waist  width’  configuration  offers  the  unique 
property  that  the  temporal  pulse  form  remains  unchanged  upon  propagation. 
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Relative  timing  stabilization.  One  part  of  the  combined  beam  is  directed  to  a  BCC 
(Supplementary  Fig.  S2),  which  consists  of  two  nearly  identical  cross-correlators 
using  200-p.m-thick  BBO  crystals,  phase-matched  for  sum-frequency  generation  of 
870  nm  light  and  2.15  pan  light.  Use  of  the  SWIR  OPCPA  delay  stage  and  a  4-mm- 
thick  calcium  fluoride  (CaF2)  window  between  cross-correlators  sets  the  group  delay 
between  pulses  to  +25  fs  in  one  cross-correlator  and  —25  fs  in  the  other.  An 
additional  2-mm-thick  calcium  fluoride  window  ensures  zero  group  delay 
(A t  =  0.0  fs)  at  the  combined  output.  For  deviations  from  this  zero-delay 
configuration  of  up  to  +20  fs,  the  photodetector  signal  is  linearly  proportional  to  the 
time  difference  and  thus  can  be  used  as  the  error  signal  fed  to  the  loop  filter  in  the 
feedback  system.  Furthermore,  in  the  vicinity  of  the  zero  crossing,  the  setup  delivers 
a  balanced  signal  and  thus  the  amplitude  noise  of  each  OPCPA  output  does  not 
affect  the  detected  error  signal. 

2DSI.  In  the  2DSI  setup  (Supplementary  Fig.  S3),  the  combined  beam  is  first  split  by 
a  beam  sampler  in  which  the  second  surface  is  anti-reflection  coated.  A  copy  of  the 
beam  (4%)  is  Fresnel-reflected  and  only  guided  via  silver  mirrors  before  being  mixed 
in  a  40  pm  type  II  BBO.  The  other  copy  of  the  beam  (96%)  passes  through  the 
beam  sampler  and  is  highly  stretched  before  being  equally  split  again  by  a  cube 
beamsplitter,  routed  to  the  BBO  and  mixed  with  the  unchirped  pulse.  Two  collinear, 
temporally  overlapped,  but  spectrally  sheared  up-converted  pulses  are  then 
generated.  To  observe  the  interference  between  the  two  up-converted  pulses,  which 
encodes  the  spectral  group  delay  information,  the  delay  of  one  of  the  highly  chirped 
pulses  is  scanned  over  a  few  optical  cycles.  The  spectrum  of  the  up-converted  signal 
is  recorded  as  a  function  of  this  delay,  yielding  a  two-dimensional  intensity  function 
that  is  shown  in  Fig.  2b, c.  The  interpretation  of  the  2DSI  data  is  relatively 
straightforward;  each  spectral  component  is  vertically  shifted  in  proportion  to 
its  group  delay. 

It  should  be  noted  that  we  treat  the  combined  beam  as  a  single  pulse,  and  use  the 
2DSI  to  retrieve  the  frequency-dependent  group  delay  of  the  synthesized  pulse,  and 
not  just  those  of  the  individual  OPCPA  pulses.  That  is,  we  measured  the  combined 
beam,  not  the  two  OPCPA  pulses  independently,  and  the  portion  mixed  with  the 
unchirped  pulse  is  purely  derived  from  the  NIR  OPCPA  such  that  the  measured 
spectral  group  delay  has  a  definite  reference  throughout  the  whole  spectrum  from 
700  to  2,500  nm.  A  different  relative  timing  results  in  a  vertical  shift  of  the  fringe 
patterns  in  Fig.  2b, c. 
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Abstract:  Noise  evolution  in  an  optical  parametric  chirped-pulse  amplifier 
(OPCPA)  differs  essentially  from  that  of  an  optical  parametric  or  a 
conventional  laser  amplifier,  in  that  an  incoherent  pedestal  is  produced  by 
superfluorescence  that  can  overwhelm  the  signal  under  strong  saturation. 
Using  a  model  for  the  nonlinear  dynamics  consistent  with  quantum 
mechanics,  we  numerically  study  the  evolution  of  excess  noise  in  an 
OPCPA.  The  observed  dynamics  explain  the  macroscopic  characteristics 
seen  previously  in  experiments  in  the  practically  important  saturation 
regime. 
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1.  Introduction 

The  optical  parametric  chirped  pulse  amplifier  (OPCPA)  [1]  is  a  ground-breaking  tool  for 
intense  laser  physics.  Combining  the  large  gain,  broad  bandwidth,  and  high  average  power 
handling  of  an  optical  parametric  amplifier  (OP  A)  with  the  high -peak-power  capability  of  a 
chirped-pulse  amplifier  (CPA),  it  is  currently  the  most  promising  technology  for  scaling  the 
peak  power  of  ultrashort  light  pulses  throughout  the  visible  [2—4]  and  near-to-mid  infrared 
spectral  ranges  [5-7].  The  technology  is  essential  both  for  high-energy,  petawatt  laser 
facilities  [8]  used  to  study  laser  fusion  [9],  proton  beam  emission  from  solids  [10],  and 
compact  optical  electron  acceleration  [11,12],  and  for  sources  of  multi-terawatt  carrier- 
envelope  phase  stable  few-cycle  pulses  used  to  control  electron  wavepacket  evolution  in 
atoms  and  molecules  [13]. 

Like  any  phase  insensitive  optical  amplifier,  OPCPAs  are  subject  to  amplified 
spontaneous  emission,  often  called  parametric  superfluorescence  (PSF),  i.e.,  parametric 
amplification  of  the  quantum  noise  due  to  two -photon  emission  from  a  virtual  level  excited  by 
the  pump  field  and  stimulated  by  the  signal  and  idler  field  zero -point  fluctuations  [14].  This 
process,  together  with  the  amplified  stimulated  emission  (ASE)  contributions  or  other 
intensity  fluctuations  of  the  pump  laser  [15,16],  affects  the  amplified  pulse  contrast.  As  we 
find  in  this  study,  the  amplification  properties  of  an  OPCPA  uniquely  affect  its  noise 
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performance.  In  a  conventional  CPA,  in  which  a  laser  medium  amplifies  a  chirped  signal 
pulse,  a  population  inversion  of  localized  emitters  provides  homogeneously  broadened  gain 
that  saturates  with  the  pulse  fluence.  In  contrast,  an  OPCPA  is  based  on  instantaneous  second- 
order  nonlinear  processes,  and  the  excited  virtual  levels  travel  with  the  pump  pulse  at  a  group 
velocity  closely  matching  that  of  the  signal;  the  virtual  levels  therefore  saturate 
instantaneously  on  the  signal’s  retarded  time  frame.  Compared  to  a  standard  OP  A,  an  OPCPA 
adds  the  complexity  of  a  map  of  instantaneous  frequency  to  temporal  coordinate  during 
amplification  due  to  the  strong  linear  chirp  which  allows  for  inhomogeneous  saturation  of  the 
available  parametric  gain. 

Experiments  have  shown  that  amplification  in  the  presence  of  PSF  in  an  OPCPA  results  in 
an  amplified  signal  field  with  two  macroscopic  components  showing  different  energy 
localization:  (i)  a  coherent  pulse  with  well  defined  temporal  chirp  matching  that  of  the 
injected  signal  pulse  (the  “seed”)  and  which  therefore  can  be  compressed  to  generate  a  Fourier 
transform-limited  pulse;  (ii)  an  incoherent  pedestal  with  phase  statistics  similar  to  that  of  PSF, 
which  remains  at  picosecond  duration  when  the  signal  pulse  is  recompressed  [17,18]. 
Henceforth,  we  refer  to  these  phenomenological  components  observed  at  the  output  of  an 
OPCPA  as  the  coherent  pulse  and  incoherent  pedestal,  respectively.  Especially  in  the  case  of  a 
broadband  signal  and  high  gain,  a  severe  impact  of  PSF  on  noise  performance  is  also  well 
documented  [5,6,19]:  it  both  degrades  the  signal  stability  and  places  an  upper  limit  on  the 
extractable  signal  energy,  due  to  transfer  of  pump  energy  to  the  incoherent  pedestal.  The 
dynamics  of  this  energy  transfer  during  amplification  have  not  been  observed  yet:  their 
understanding  in  the  highly  nonlinear  saturation  regime  is  not  only  of  fundamental  interest, 
but  is  particularly  important  with  regard  to  performance,  since  amplifier  saturation  is 
necessary  for  obtaining  good  conversion  efficiency  and  stable  output  energy.  Intuition,  based 
on  the  properties  of  laser  and  electrical  amplifiers,  suggests  that  saturation  should  suppress 
fluctuations.  However,  the  effects  of  saturation  on  excess  quantum  noise  lack  a  satisfactory 
description:  while  the  experiments  show  degradation  of  pulse  contrast  during  saturation 
[6,19],  a  numerical  analysis  of  output  statistics  of  PSF  in  an  OPCPA  seeded  by  a  distributed 
classical  noise  source  did  not  isolate  the  effect  [20]. 

In  this  paper,  we  introduce  a  quantum-mechanically  consistent  numerical  model  of  the 
dynamics  of  PSF  growth  in  an  OPCPA  that  captures  the  process  of  energy  exchange  during 
amplification  between  what  will  become  coherent  and  incoherent  components  of  the  electric 
field  after  compression,  and  well  reproduces  the  macroscopic  characteristics  observed  in 
experiments.  Since  the  purpose  of  this  paper  is  to  isolate  the  influence  of  PSF  on  the  amplified 
pulse  contrast,  we  do  not  consider  effects  due  to  ASE  contributions  or  other  intensity 
fluctuations  of  the  pump  laser.  By  virtue  of  the  model’s  adherence  to  quantum  mechanics 
even  in  the  highly  nonlinear  saturation  stage,  we  observe  the  saturation  dynamics  of  a 
quantum-noise-contaminated  OPCPA,  uncovering  several  distinguishing  features.  We  find 
that  PSF  must  be  characterized  by  two  observables  which  display  different  evolution 
dynamics:  the  shot-to-shot  energy  fluctuation  and  the  ratio  of  coherent  pulse  energy  and 
incoherent  pedestal  energy  of  the  amplified  signal  field.  We  find  that  an  OPCPA  has  well 
defined  but  different  operating  points  for  maximum  suppression  of  PSF-induced  fluctuations 
or  pedestals.  Beyond  these  operating  points,  heavy  saturation  leads  to  large  excess  noise  that 
can  be  enhanced  by  orders  of  magnitude. 

2.  Numerical  model 

For  the  numerical  description  of  the  PSF  dynamics  in  the  amplification  process,  we  focus  on 
an  OPCPA  seeded  by  the  initial  quantum  noise  field  and  a  chirped  signal  field;  amplification 
occurs  in  a  periodically-poled  stoichiometric  lithium  tantalate  (PPSFT)  crystal.  The  nonlinear 
quantum  system  dynamics  can  be  described  by  a  quasi-probability  distribution,  such  as  the 
Wigner  distribution  (WD)  [21].  It  is  well  known  from  quantum  optics  that  for  linear  systems 
the  evolution  equation  for  the  WD  is  equivalent  to  a  classical  Fokker-Planck  equation,  and  is 
thus  also  equivalent  to  a  stochastic  process  involving  classical  noise  sources,  resulting  in  a 
semiclassical  picture  of  the  quantum  process.  This  correspondence  has  been  exploited  in 
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numerical  studies  of  PSF,  OPA,  and  optical  parametric  oscillation  in  their  linear  regimes 
[22,23].  It  is  less  well  known  that  for  the  case  of  weak  nonlinearities,  i.e.,  no  significant 
nonlinear  effects  at  the  few-photon  level,  the  Fokker-Planck  approximation  holds  and  the 
nonlinear  quantum  system  dynamics  can  still  be  extracted  accurately  from  stochastic 
Langevin  equations  [24,25],  an  approach  used  earlier  to  study  the  quantum  noise  in  parametric 
amplifiers  used  for  squeezed  light  generation  [26,27].  These  stochastic  equations  have  a 
deterministic  component  equal  to  the  Heisenberg  equations  of  motion  for  the  field  operators 
and  are  complemented  by  relaxation  terms  and  associated  noise  terms.  For  the  case  of  a 
lossless  OPA  process,  fluctuations  stem  solely  from  the  quantum  mechanical  uncertainty  in 
the  input  fields.  Knowledge  of  a  quasi-probability  distribution  allows  computation  of  all 
expectation  values  of  quantum  mechanical  observables,  and  for  the  case  of  the  WD  and  its 
associated  stochastic  process,  computed  expectation  values  correspond  to  quantum 
mechanical  expectation  values  of  symmetrically  ordered  field  operators  [21].  Thus,  this 
approach  allows  for  a  rigorous  treatment  of  quantum  fluctuations  in  weakly  nonlinear 
quantum  optical  systems,  such  as  OPCPAs  with  large  mode  cross  sections.  For  completeness, 
we  also  mention  that  the  quantum  dynamics  of  a  second-order  nonlinear  process,  as  is  the  case 
discussed  here,  can  be  described  exactly  with  the  help  of  the  positive  P-representation, 
pioneered  by  Peter  Drummond  [28].  It  has  been  shown  [29]  that  third-  and  higher-order 
moments  of  the  electric  field  can  differ  significantly  whether  calculated  by  means  of  the  P  or 
truncated  WD  representation.  However,  for  large  normalized  photon  numbers  and  weak 
nonlinearity,  differences  are  very  small  (the  discrepancy  was  quantitatively  small  for  pump 
photon  numbers  of  100  in  [25],  though  some  slight  differences  persisted  at  early  times.).  With 
the  gigawatt  peak  powers  typical  of  OP  As,  we  are  always  in  the  high  photon  number  (totaling 
~1015  in  our  case)  and  weak  nonlinearity  limit.  This  fact,  together  with  the  increased 
mathematical  complexity  of  the  positive  P-representation,  using  a  twice  as  large  phase  space 
which  considerably  increases  computation  time,  led  us  to  work  within  the  truncated  WD. 

We  simulate  the  evolution  of  noise  in  an  OPCPA  by  numerically  solving  the  coupled 
nonlinear  equations  of  parametric  amplification  in  the  spectral  domain,  accounting  for  linear 
dispersion  to  all  orders  [30,31];  the  equations  describe  the  interaction  among  signal,  idler  and 
pump,  which  in  the  following  will  be  respectively  labelled  i  =  1,  2,  3.  These  waves  propagate 
along  the  z  coordinate  with  carrier  frequency  ojt  and  wavenumber  kt.  To  exploit  the  large  d33 
nonlinear  coefficient  of  the  crystal  for  all  fields  polarized  along  the  extraordinary  axis,  we 
operate  in  the  quasi  phase-matching  regime,  obtained  by  periodically  poling  the  nonlinear 
crystal:  poling  is  accounted  for  by  changing  the  sign  of  along  z-  The  carrier  fields  therefore 
experience  at  any  crystal  coordinate  a  real  phase  mismatch  A k  =  krk2-kh  We  describe  the 
electric  field  of  each  wave  as: 

E, (z,t)  =  U2-{Ai(z,t) exp [j(oj2 - ktz)]  +  c.c.}  =  sl?{A(z,f) exp [jio^t - kiz)])  (1) 


where  Afat)  denotes  the  field  complex  amplitude.  The  coupled  equations  describing  the 
second  order  interaction  of  the  fields  are  derived  from  the  nonlinear  propagation  equation: 


d2E  s2dl  d2pNL 
W-Po^r  =  M0^r- 


(2) 


applied  on  the  total  field  E(z,t )  =  Ej(z,t)  +  E2(z,t )  +  E3{z,t). 

Here  DL(z,t )  =  s0^sr(T)E(z,t-T)dr  is  the  linear  electric  induction  field  accounting  for  the 

linear  dispersion  of  the  medium  [32],  and  PNL  =  2s0deff  E(z,t)2  is  the  nonlinear  polarization, 
where  de ff  is  the  effective  second-order  non-linear  coefficient.  Since  our  model  also  accounts 
for  a  broadband  noise  field,  and  our  purpose  is  to  calculate  the  evolution  of  the  noise  with  the 
highest  accuracy,  we  avoid  the  slowly-varying-envelope  approximation  [33]  typically  adopted 
to  simplify  the  calculations;  in  addition  we  consider  the  linear  dispersion  of  the  material  to  all 
orders.  For  this  purpose,  we  develop  Eq.  (2)  in  the  frequency  domain  by  taking  its  Fourier 
transform  and  obtaining: 
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(3) 


d2E  Q2n(Q)2 


E  —  /j0Q.  PN1 


where  E(z,Q)  =  E{E(z,t)}  and  PM(z,Cl)  =  e{Pn,(zJ))  are  the  non-unitary  Fourier 

transforms  of  the  electric  field  and  of  the  nonlinear  polarization  respectively,  Q  is  the  angular 
frequency  and  n(Q)  is  the  frequency-dependent  refractive  index.  PNL  is  developed  rejecting 
components  at  frequencies  different  from  ©i,  co2  and  co3;  when  the  fields  Et  are  not  overlapped 
in  frequency,  wave  vector,  and  polarization,  it  is  possible  to  split  Eq.  (3)  into  three  coupled 
equations  which  separately  describe  the  evolution  of  the  field  envelopes: 


9  A  d\  ,27  -jAkz 

—j--j2kl—±  +  blAl=-cl-e 

OZ  oz 

3^2  *0  7  7  2  7  _  -jAkz 

~  2  2  A  A  —  C2  ’  e 

dz  oz 

*A-j2k,d-^  +  brA,  =  -c,-e+>M’- 


(4) 


dz 2  '  "  dz 

In  this  case  Ai{z,co)  =  (F{Ai(z,t )}  is  the  Fourier  transform  of  the  envelope  amplitude  of 

each  field  and  co  =  Q  -  is  the  detuning  from  the  carrier  frequency  oj,\  coefficients  l?~  and  c, 
are  defined  as: 


b2  = 


-k2  +ki 2 ,  with  Z,  -  +  •»,-(<»  + <**) 


(5a) 


and 


f  ^2 

CO  +  cox 


CQ+CQ- 


dtff ^{AA  *}  ’ 


(5b) 


Co  = 


y 

'3 


Aff^jAA} 


V  "0  J 

Here  n^co  +  co?)  are  the  refractive  index  functions  deduced  from  the  Sellmeier  equations, 
and  allow  to  take  into  account  the  whole  linear  dispersion  of  the  material.  The  system  can  be 
solved  as  follows:  if  a  suitably  small  step  Az  is  chosen,  the  products  A2A2*,  A2A7*  and  A7A2 
are  nearly  constant,  and  Eqs.  (4)  can  be  analytically  solved.  Given  the  fields  ^{z.co)  at  the 
beginning  of  a  step,  the  fields  at  z  +  Az  are: 


A(z  +  Az,o)  * 

Aj(z,<w)  h — - 

•  cxpl  /T/c,  -  Zj)AzJ—  —  •  expf-  jAkAz] 

n_ 

Yi 

A,(z  +  Az,a>)  * 

A,(z,®)  +  — 
Yi. 

■  exp[/(fc2  -  £2)Azl-  —  •  expf-  JAkAz] 

Y 2 

Aj(z  +  Az,co)  a 

A3(z,co)  +  — 

L  r3_ 

•  expf j(k3  -  k3)Az  -  —  •  exp[+  jAkAz] 

Yi 

where  yl2  =bl2  -Ak2  -2 kl2Ak  and  y3  =  b3  -Ak2+2k3Ak  .  It  is  important  to  remark  that 
the  representation  of  fields  given  in  Eqs.  (4)  assumes  that  pump,  signal  and  idler  are  three 
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separate  fields.  Justifying  this  treatment,  the  amplifier  we  model  employs  a  small  non- 
collinear  angle  between  signal  and  idler,  used  both  to  allow  their  separation  after 
amplification  (since  they  have  opposite  temporal  chirp)  and  to  avoid  signal-idler  interference 
for  preservation  of  carrier-envelope  phase  of  the  signal.  The  results  we  obtain  from  our  model 
do  not  hold  for  degenerate  collinear  OP  As,  for  which  fields  1  and  2  of  Eqs.  (4)-(6)  collapse 
into  one  equation,  and  for  strongly  non-collinear  geometries,  which  would  require  at  least  one 
more  spatial  coordinate. 

Our  1-D  plane  wave  model  includes  all  longitudinal  modes,  m,  and  their  associated  noise. 
In  the  frequency  domain,  at  any  mode  frequency,  com,  the  corresponding  component  of  the 
initial  signal,  idler,  or  pump  electric  field  is  represented  by  a  complex  stochastic  phasor,  Am(0) 
=  Bm  +  nm.  Bm  is  the  deterministic  component  of  the  field,  and  is  set  to  0  in  the  case  of  the 
idler;  nm  is  a  zero-mean,  stochastic  phasor  representing  the  independent  fluctuations  of  the 
field. 


Real  part 


Fig.  1.  (Color  online)  (a)  Schematic  representation  of  the  deterministic  part  (arrow)  and  50 
stochastic  components  (circles)  of  a  field  mode,  (b)  Initial  signal  field  distribution  (scatter), 
evaluated  at  three  modes  of  frequency  com,  experiencing  respectively  the  highest  gain  G0,  G0/2 
and  Go/10  and  separated  by  a  phase  shift  imparted  by  the  chirp,  (c)  Depiction  of  the  same 
modes  after  amplification  and  (d)  after  compression;  pedestal  fields  deduced  after  subtracting 
the  deterministic  components  are  indicated  in  the  dashed  circle.  All  data  refer  to  configuration 
II.  Panels  (a)  and  (b)  are  magnified  500  times  with  respect  to  (c)  and  (d). 

Note  that  fluctuations  are  included  for  each  of  the  signal  and  idler  fields;  the  quantum 
noise  of  the  pump  is  negligible,  as  we  also  confirmed  independently  by  simulations  not 
reported  here.  Real  and  imaginary  components  of  nm  are  taken  as  uncorrelated  Gaussian 
distributions  [22]  with  variance  am2  oc  com.  A  representation  of  these  fluctuating  fields  is  given 
in  Figs.  1(a)  and  1(b),  where  we  show  Bm  (vectors),  and  Am  and  nm  (scatter)  for  three  modes 
com  of  the  signal  field.  Our  numerical  method  treats  identically  the  initial  electric  field 
components  whether  originating  from  the  deterministic  field  or  the  vacuum  fluctuations.  We 
apply  this  tool  to  the  study  of  an  ultra-broadband  OPCPA  system  known  to  be  sensitive  to 
PSF  [5,6].  We  model  a  typical  high-gain  pre-amplifier,  in  which  noise  begins  at  the  vacuum 
level  and  that  establishes  the  noise  content  of  later  stages  of  a  multi-stage  system.  The 
amplifier,  pumped  by  a  9-ps  FWHM  Gaussian  pulse  at  1.047  pm  and  seeded  by  a  broadband 
(69-THz  FWHM)  pulse  at  2.094  pm  for  operation  around  degeneracy,  uses  a  3 -mm  long 
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PPSLT  crystal  with  poling  period  A  =  31.2pm.  These  parameters  are  close  to  the  experimental 
conditions  of  Ref.  [6].  The  pump-to-seed  energy  ratio  is  106. 


Table  1.  Key  Parameter  for  Three  OPCPA  Configurations" 


Input  (z 

=  0  mm) 

Output  (z  =  3  mm) 

a 

o 

Seed  duration 
[ps] 

Pump 

intensity 

[GW/cm2] 

n  [%] 

Av 

[THz] 

rpAv 

[%-THz] 

SPR 

SNR 

I 

1.04 

6.254 

9.3 

56.8 

528.2 

17.1 

310.6 

II 

7.34 

8.356 

48.5 

43.5 

2111.9 

54.2 

412.7 

III 

10.48 

8.825 

56.6 

27 

1528.2 

52.6 

322.2 

"For  efficiency  r|,  bandwidth  Av  and  efficiency-bandwidth  product  columns,  maximum  values  are  highlighted  in  bold. 

The  variances  am2  are  determined  by  the  quantum  fluctuations  due  to  the  longitudinal 
modes  of  the  100-ps-long  simulation  window.  This  number  is  further  increased  by  a  factor 
equal  to  the  number  of  transverse  modes  amplified  assuming  a  pump  beam  of  100-pm  radius 
in  the  3 -mm  long  crystal,  which  is  estimated  as  about  25.  We  note  that  this  choice  results  in  a 
calculated  amplified  pulse  contrast  that  closely  matches  that  observed  in  equivalent 
experiments  [6].  This  system  was  investigated  previously  in  order  to  find  the  conditions  of 
maximum  efficiency-bandwidth  product  [34].  Following  that  analysis,  we  choose  3  values  of 
seed  chirp,  summarized  in  Table  1,  which  correspond  to:  an  under-chirped  amplifier,  with 
maximum  amplified  signal  bandwidth  but  limited  conversion  efficiency  (Configuration  I);  an 
amplifier  chirped  for  maximum  efficiency-bandwidth  product  (Configuration  II);  and  an  over¬ 
chirped  amplifier  (Configuration  III),  with  excellent  conversion  efficiency  but  significant 
spectral  narrowing.  Table  1  also  provides  the  pump  intensity  corresponding  to  each  of  the 
three  regimes.  For  each  configuration,  we  evaluated  50  independent  trajectories  triggered  by 
uncorrelated  noise  fields.  The  averages  taken  over  this  ensemble  of  classical  solutions 
correspond  to  quantum-mechanical  expectation  values  [35].  The  results  of  a  batch  of 
simulations  are  depicted  in  Fig.  1(c). 
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Frequency  [THz] 


Fig.  2.  (Color  online)  WD  map  of  a  signal  field  from  configuration  II,  evaluated  before  (a)  and 
after  (b)  compression. 

The  incoherent  nature  of  the  amplified  noise  is  evidenced  in  panel  (d),  representing  the 
same  fields  after  compression.  A  complete  description  of  an  amplified  signal  field  is  given  in 
Fig.  2,  showing  the  WD  before  (a)  and  after  (b)  compression.  The  WDs  clearly  reveal  the 
presence  of  a  strong  incoherent  field,  arising  from  PSF,  superimposed  to  the  coherent  chirped 
amplified  signal;  as  a  feature  common  to  all  configurations,  this  component  has  a  spectrum 
corresponding  to  the  phase -matching  bandwidth  of  the  OPCPA,  and  duration  comparable  to 


#142105  -  $15.00  USD  Received  8  Feb  2011;  revised  30  Mar  2011;  accepted  30  Mar  2011;  published  15  Apr  2011 
(C)  201 1  OSA  25  April  201 1  /  Vol.  19,  No.  9  /  OPTICS  EXPRESS  8363 


the  pump  pulse.  After  applying  a  linear  chirp  to  compress  the  coherent  pulse  (see  Fig.  2(b)), 
the  incompressible  temporal  pedestal  degrades  the  pulse  contrast,  as  observed  in  experiments. 

In  addition  to  this  pedestal,  PSF  also  strongly  affects  the  shot-to-shot  energy  stability  of 
the  amplified  signal.  The  noise  contaminating  the  signal  field  during  the  OPCPA  process  can 
thus  be  characterized  by  two  quantities:  (i)  the  signal -to-noise  ratio  (SNR),  which  evaluates 
the  shot-to-shot  energy  fluctuation  of  the  signal  pulse;  (ii)  the  signal -to-pedestal  ratio  (SPR), 
which  measures  the  ratio  between  the  energies  of  the  coherent  pulse  and  the  incoherent 
pedestal.  To  calculate  the  SNR  we  must  derive  the  expectation  value  and  variance  of  a  mode 
intensity  Im(z)  =  amf(z)  am(z)  of  signal  or  idler,  (where  am|(z),  am(z)  are  the  photon  creation 
and  annihilation  operators),  in  terms  of  expectation  values  of  the  classical  stochastic  field 
variables  Am(z),  Am*(z).  Equating  the  computed  expectation  values  of  the  stochastic  field 
variables  to  the  expectation  values  of  the  corresponding  symmetrically  ordered  quantum 
mechanical  field  operators,  we  derive: 

(4U))  =  (k(2)f)-(k|2)  (7a) 

kmCz)!  =  (Akoo \2J  -  (kf)  (7b) 

Equations  (7)  give  us  a  physically  satisfying  answer  for  the  initial  field  intensity  and 
variance:  (/m(0))  =  |i?m|2 ,  /.<?.,  the  deterministic  component  of  the  initial  field,  and 

[A/m(0)]2  =  2\Bmf  ,  z.e.,  the  corresponding  photon  number  fluctuations  or  shot  noise. 

These  equations  also  correctly  give  a  null  expectation  value  and  null  variance  for  a  vacuum 
state  (Bm  =  0):  this  is  of  key  importance,  since  the  simulation  is  initially  dominated  by  modes 
in  the  vacuum  state.  From  these  quantities  we  can  calculate  the  expectation  value  (E(z))  and 

variance  A E(z)  of  the  pulse  energy: 


E(z)=^M)Aa,m  (8a) 

m 

[A E{zf=((E(z)-{E(z))f)  (8b) 

and  define  SNR(z )  =  (£’(z))/AE(z)  •  Figure  3  shows  the  evolution  with  propagation  length  of 

(E(z))  and  A E(z)  for  signal  and  idler,  while  Fig.  4(a)  shows  SNR(z)  (solid  lines).  The 

pedestal  field  is  deduced  from  the  comparison  between  the  PSF -contaminated  signal  field 
Am(z)  and  a  reference  field  ARm(z )  obtained  from  the  simulation  of  a  noise-free  amplifier  (nm  = 
0).  In  Figs.  1(c)  and  1(d)  the  dots  correspond  to  the  PSF-contaminated  signal  fields  and  the 
arrows  to  the  deterministic  reference  fields.  We  find  that  PSF  does  not  significantly  modify 
the  amplification  regime,  in  terms  of  pump  depletion;  the  pedestal  field  APm(z )  can  be 
therefore  evaluated  as  APm(z)  =  Am(z)-ARm(z)-  The  SPR  is  finally  evaluated  by  taking  the  ratio 
between  the  signal  and  the  pedestal  pulse  energy. 

It  should  be  noted  that,  while  the  SNR  can  be  measured  by  the  shot-to-shot  fluctuations  of 
the  amplified  pulses,  the  SPR  is  difficult  to  characterize  experimentally.  In  fact,  the  effect  of 
PSF  cannot  be  isolated  by  measuring  an  unseeded  OPCPA,  since  pump  depletion  is  much 
lower  than  in  the  seeded  case.  The  experimental  measurement  of  the  evolution  of  these 
quantities  during  the  OPCPA  process  poses  an  even  greater  challenge.  Our  numerical 
approach,  on  the  other  hand,  has  the  unique  capability  of  isolating  the  pedestal  field  from  the 
signal  under  realistic  conditions,  allowing  exploration  of  its  dependence  on  amplifier 
parameters,  as  well  as  SNR  and  SPR  evolution  in  the  OPCPA  crystal. 
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3.  Results  and  discussion 


In  Fig.  3  we  report  the  evolution  of  energy  mean  and  standard  deviation  of  signal  and  idler, 
for  configurations  I  and  II;  for  completeness,  we  have  extended  the  calculation  beyond  the 
optimum  crystal  length  of  3  mm.  From  the  curves  we  can  clearly  identify  three  stages  of 
amplification:  growth,  saturation,  and  over-saturation. 
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Fig.  3.  (Color  online)  Evolution  of  energy  mean  and  standard  deviation  for  signal  and  idler;  the 
energy  growth  of  the  pedestal  mean  is  also  given.  In  the  case  of  configuration  II  we  indicate  the 
coordinate  at  which  the  pump  peak  is  fully  depleted  (see  inset).  Trends  calculated  for 
configuration  III  (not  shown)  are  comparable  to  the  ones  of  configuration  II. 

At  the  start  of  amplification,  the  unseeded  idler  intensity  catches  up  with  the  signal  and  the 
signal  standard  deviation  grows  relative  to  the  mean  as  a  result  of  mixing  with  the  idler, 
causing  a  degradation  of  the  initial  SNR  [Fig.  4(a)].  After  this,  exponential  growth  of  both 
waves  sets  in.  In  this  regime  we  observe  a  gradual  decay  of  SNR. 
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Fig.  4.  (Color  online)  (a)  Evolution  of  SNR  and  SPR  for  the  three  configurations  of  Table  1 .  (b) 

Temporal  profiles  of  the  signal  pulses  after  compression,  each  normalized  to  its  peak. 

Once  we  enter  the  pump  depletion  regime,  however,  the  signal  energy  fluctuation  is 
dramatically  reduced,  thus  improving  the  SNR.  The  SNR  reaches  a  sharp  maximum  shortly 
before  the  exit  facet  of  the  crystal,  close  to  the  propagation  distance  for  which  the  pump  peak 
fully  depletes  and  the  maximum  pump  to  signal/idler  conversion  occurs.  After  this  point,  a 
rapid  drop  in  SNR  is  seen.  This  clearly  shows  that,  while  in  the  pump  depletion  regime  the 
SNR  can  be  considerably  enhanced,  in  over- saturation  the  mean  intensity  sags  due  to  back- 
conversion  to  the  pump,  while  noise  continues  to  grow  on  average  across  the  pulse.  Figure 
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4(a)  also  shows  that,  while  the  three  configurations  exhibit  very  similar  trends  up  to 
saturation,  configuration  II  is  preferable  due  to  its  higher  SNR  value. 

After  saturation,  the  low-chirp  configuration  results  in  the  most  dramatic  drop  in  SNR,  a 
result  of  strong  growth  of  PSF  in  the  unseeded  temporal  wings  while  the  amplified  signal 
converts  back  to  pump  at  the  pulse  center.  In  comparison,  the  SPR  trend  reveals  dramatic 
differences  with  that  of  SNR,  indicating  that  the  energy  of  the  incoherent  pedestal,  evaluated 
through  the  SPR,  is  not  directly  correlated  to  the  pulse  energy  fluctuations  measured  by  the 
SNR.  In  fact,  maximum  suppression  of  energy  fluctuations  in  configuration  I  occurs 
simultaneously  with  strong  degradation  of  pulse  contrast,  and  in  this  case  the  SPR  is  sharply 
reduced  at  the  peak  of  saturation,  i.e.,  at  z  =  3  mm,  where  conversion  efficiency  and 
efficiency-bandwidth  product  are  maximized.  The  two  parameters  therefore  offer 
complementary  information  for  evaluating  the  impact  of  excess  noise  in  the  OPCPA.  The 
profile  of  the  incoherent  pedestal  after  pulse  compression  is  given  in  Fig.  4(b);  as  expected, 
configuration  I  is  strongly  contaminated  by  an  incoherent  pedestal  generated  at  unseeded 
temporal  coordinates.  Note  that  the  different  pedestal  durations  are  due  to  pulse  compression, 
which  imparts  a  vertical  shearing  to  the  WDs  (see  Fig.  2).  These  observed  dynamics  explain 
the  experimental  results  of  Ref.  [6],  including  the  non-intuitive  result  that  a  saturated  amplifier 
can  exhibit  reduced  pulse  contrast.  We  note,  the  observed  degradation  of  SPR  in  saturation, 
and  of  both  SPR  and  SNR  in  oversaturation,  are  characteristic  phenomena  singular  to 
quantum  noise  growth  in  an  OPCPA.  In  a  conventional  CPA,  there  is  no  oversaturation 
regime;  while  nonlinear  wave  mixing  in  an  OPCPA  allows  back-conversion  of  signal  and 
idler  to  pump  after  amplification,  saturation  of  a  population  inversion  in  a  laser  is  not 
reversible.  Laser  gain  saturation  is  also  homogeneous,  and,  therefore,  unlike  the  OPCPA, 
unseeded  frequency  modes  of  the  signal  field  cannot  continue  to  amplify  exponentially  while 
the  gain  for  seeded  modes  saturates.  In  a  conventional  OP  A,  the  absence  of  a  chirp  ensures 
that  saturation  and  back-conversion  is  uniform  with  respect  to  frequency,  and  no  macroscopic 
pedestal  is  produced,  thus  allowing  SNR  to  fully  characterize  the  noise  performance. 

4.  Conclusions 

In  conclusion,  we  have  performed  a  quantum-mechanically  consistent  numerical  investigation 
of  the  dynamics  of  PSF  growth  in  a  realistic  high-gain  OPCPA.  Thanks  to  the  model’s 
capability  to  simulate  also  the  saturation  stage  of  amplification,  this  investigation  for  the  first 
time  captures  all  dynamics  of  a  quantum-noise-contaminated  OPCPA,  addressing  the 
important  practical  issues  of  signal  energy  stability  and  pulse  contrast.  Both  quantities  are 
related  to  the  incoherent  pedestal  superposed  to  the  signal,  but  display  different  evolution 
dynamics  throughout  the  amplification  process.  Three  operating  conditions  were  explored, 
characterized  by  different  chirps  of  the  input  seed  and  maximizing,  respectively,  the 
bandwidth,  the  efficiency-band  width  product,  and  the  conversion  efficiency.  We  find  that  the 
chirp  maximizing  the  efficiency-bandwidth  product  is  also  characterized  by  the  smallest 
contribution  of  the  noise,  both  in  terms  of  energy  fluctuation  and  of  pulse  contrast. 
Significantly,  we  find  that  while  amplifier  saturation  improves  the  signal’s  shot-to-shot  energy 
stability,  it  does  not  necessarily  improve  the  pulse  contrast.  In  fact,  for  the  case  of  an  amplifier 
optimized  for  bandwidth,  strong  degradation  of  the  pulse  contrast  (i.e.,  growth  of  the  PSF- 
derived  incoherent  pedestal  energy  relative  to  the  coherent  signal  pulse  energy)  is  observed  in 
saturation.  Over- saturation,  sometimes  employed  to  boost  amplifier  bandwidth  at  the  expense 
of  conversion  efficiency,  uniformly  reduces  both  pulse  contrast  and  energy  stability. 

Knowledge  of  these  dynamics  thus  provides  important  insight  for  the  optimization  of 
OPCPA  systems  applied  to  the  study  of  strong-field  laser  physics,  as  well  as  increases  our 
fundamental  understanding  of  quantum  noise  in  parametric  amplification. 
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Abstract:  We  present  a  scheme  for  the  generation  of  ultrabroadband  light  pulses  based  on  the  coherent  synthesis  of  the 
pulses  delivered  by  two  broadband  optical  parametric  amplifiers.  The  scheme  provides  pJ-level  pulses  with  a  spectrum 
ranging  from  500  to  1000  nm,  supporting  sub-4  fs  duration.  The  preliminary  results  obtained  here  demonstrate  that  this 
scheme  is  suitable  for  the  generation  of  single-cycle  CEP-stable  pulses,  whose  energy  can  be  boosted  by  additional 
amplification  stages. 


Few-cycle  optical  pulses  are  currently  used  to  drive 
strong-field  processes  such  as  high-harmonic  generation 
(HHG),  which  has  recently  led  to  the  generation  of 
isolated  attosecond  pulses  as  short  as  80-as  in  the  XUV 
range1.  These  processes  are  sensitive  to  the  electric  field 
of  the  pulse  rather  than  its  envelope  profile,  thus 
requiring  control  of  the  pulses  Carrier-Envelope  Phase 
(CEP)  to  obtain  a  reproducible  electric  field.  In 
addition,  a  broadband  source  is  necessary  to  access  the 
few-to-single  optical  cycle  regime.  White-light  seeded 
Optical  Parametric  Amplifiers  (OPAs)  are  powerful 
tools  for  the  generation  of  ultrabroadband  pulses,  since 
they  can  provide,  under  suitable  phase-matching 
conditions,  ultrabroad  gain  bandwidths2.  In  addition 
OPAs  offer  the  opportunity  to  produce  pulses  with 
passively  stabilized  CEP,  exploiting  the  difference 
frequency  generation  process  that  occurs  in  the  idler 
beam  synthesis3.  Broad  gain  bandwidth  in  an  OPA  is 
achieved  when  the  group  velocities  of  signal  and  idler 
are  matched2;  this  condition  is  satisfied  either  in  the 
case  of  type  I  phase  matching  at  degeneracy,  or  in  the 
non-collinear  OPA  (NOP A),  in  which  the  idler  group 
velocity  is  projected  along  the  signal  propagation 
direction.  Using  these  concepts,  a  variety  of  broadband 
OPA  schemes,  pumped  by  either  the  fundamental 
frequency  (FF)  or  the  second  harmonic  (SH)  of 
Ti: Sapphire,  have  been  demonstrated.  In  particular,  6-fs 
visible  pulses  can  be  routinely  generated  by  SH-pumped 
NOP  A4,  and  sub-7-fs  pulses  are  produced  in  the  650-nm 
to  950-nm  spectral  region  by  the  SH-pumped 
degenerate  OPA5.  Since  the  WLC  presents  a  highly 
structured  intensity  and  phase  profile  around  the  driving 
pulse  wavelength,  for  the  degenerate  OPA  additional 
stages  are  necessary  to  shift  the  WLC  driving 
wavelength  to  the  infrared  and  to  generate  well-behaved 
WLC  around  800  nm. 

In  this  work  we  present  a  scheme  for  the  generation 
of  CEP-stable  (sub) single-cycle  pulses  in  the  spectral 
range  from  500  to  1000  nm.  The  scheme,  summarized 


in  the  block  diagram  of  Fig.  1,  combines  (i)  the 
generation  of  a  passively  CEP-stabilized  WLC  seed 
with  (ii)  two  broadband  OPAs  seeded  by  distinct 
portions  of  the  WLC  and  (iii)  coherent  combination  of 
their  outputs  in  order  to  synthesize  high  energy 
(sub)single-cycle  pulses  with  controlled  electric  field 
profile. 

The  system  is  powered  by  550-pJ  pulses  from  an 
amplified  Ti:sapphire  laser,  providing  100-fs  pulses  at 
800-nm  wavelength  and  1kHz  repetition  rate.  The  first 
block5  is  a  two-stage  IR  OPA  for  the  generation  of  high 
energy  IR  pulses.  When  pumped  by  250-pJ  pulses,  this 
OPA  provides  a  total  signal+idler  energy  of  35  pJ. 


Fig.  1:  setup  for  separate  amplification,  compression 
and  coherent  combination  of  visible  and  near-infrared 
pulses.  SHG:  second  harmonic  generation;  WLG:  white 
light  generation. 

Both  the  signal  and  the  CEP-stable  idler  can  be 
used  to  generate  a  WLC  by  self-phase-modulation 
either  in  a  3 -mm  thick  sapphire  plate  or  in  a  4-mm  thick 
YAG  respectively.  This  broadband  seed  ranges  from 
500  to  1000  nm,  and  covers  without  gaps  the 
amplification  range  of  the  subsequent  visible  NOPA 
and  SH-pumped  degenerate  OPA  (see  Fig  2).  For  the 
preliminary  demonstration  of  the  overall  scheme 
reported  here,  we  tuned  the  OPA  at  1.3  pm  and  we  used 
its  signal  beam  for  WLC  generation  in  sapphire.  The 
supercontinuum  is  then  split  by  a  neutral-density 


ultrathin  beam-splitter,  and  the  two  replicas  are 
amplified  by  the  visible  NOPA  and  by  the  degenerate 
OPA.  The  OP  As  are  pumped  by  the  SH  from  the 
remaining  300-pJ  FF  pulses;  thick  SH  crystals  were 
used  to  narrow  the  bandwidth  of  the  SH  pulse,  and  to 
facilitate  its  temporal  overlap  with  the  chirped  WLC. 
Each  of  the  two  OPAs  provides  pulses  with  energy  of  1- 
2  pJ. 

To  ensure  good  spatial  overlap  of  the  collimated 
beams,  the  two  OPAs  are  designed  in  order  to  employ 
the  same  focal  lengths  and  to  propagate  the  beams  for 
the  same  distances.  Compression  is  separately  achieved 
by  visible  and  IR  chirped  mirrors,  which  shorten  the 
pulse  duration  close  to  the  transform-limit. 
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Fig.2:  Spectrum  of  the  gap-free  supercontinuum  (blue 
solid  line)  generated  in  the  Y AG/sapphire  plate, 
compared  to  the  gain  of  the  visible  and  infrared  OPAs 
that  will  be  employed  for  its  amplification. 


The  two  amplified  pulses  are  then  synchronized  by 
a  delay  line  equipped  with  a  piezoelectric  actuator,  and 
collinearly  combined  by  a  second  ultrathin  beam¬ 
splitter.  A  gap-free  spectrum  arising  from  the 
combination  of  the  two  pulses  is  shown  in  Fig.  3;  it 
supports  sub-4  fs  pulse  duration. 


Wavelength  (nm) 

Fig.3:  Spectrum  of  the  synthesized  pulses;  the  spectrum 
supports  sub-4  fs  pulse  duration. 

The  last  challenging  step  of  the  pulse  synthesis  is 
the  coherent  combination  of  the  two  pulses,  which  calls 
for  careful  control  of  their  relative  delay  and  phase. 
Thanks  to  the  spectral  overlap  of  the  two  beams,  we 


could  directly  use  spectral  interferometry  to  characterize 
their  delay  fluctuations;  for  this  purpose,  the  delay 
between  the  unlocked  pulses  was  increased  to  about  200 
fs  and  the  spectral  fringes  were  monitored  for  one 
minute,  giving  slow  delay  fluctuations  with  rms  of  the 
order  of  3  fs.  When  seeded  with  CEP  stable  continuua, 
we  plan  to  lock  the  relative  delay  of  the  pulses  by  using 
a  balanced  nonlinear  cross-correlator6,  which  allows 
attosecond-precision  relative  timing  stability  thanks  to 
the  capability  of  the  balanced  detection  scheme  to 
cancel  the  amplitude  noise7.  This  scheme,  together  with 
the  possibility  to  seed  the  OPAs  with  CEP-stable  light, 
will  allow  the  generation  of  single-cycle  pulses  with 
fully  controlled  electric  field.  An  additional  challenge 
will  be  full  temporal  characterization  of  the  synthesized 
ultrabroadband  pulse,  which  will  be  performed  by  two- 
dimensional  spectral-shearing  interferometry8. 
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In  the  past  several  years,  a  new  class  of  ultrashort  pulse  sources  for  long-wavelength-driven  strong-field 
experiments  has  been  developed  based  on  ultrabroadband  optical  parametric  chirped  pulse  amplification 
(OPCPA)  with  2-pm  signal  wavelength  [1-3].  These  amplifiers  share  the  advantageous  properties  of  few-cycle 
duration,  passive  carrier-envelope  phase  control  and  close-to-millijoule  energy,  but  also,  due  to  low  seed  energy, 
high  gain,  and  large  mode  cross  section,  a  propensity  towards  significant  contamination  from  amplified  vacuum 
noise.  This  amplified  superfluorescence  has  been  found  both  to  degrade  the  signal-to-noise  ratio  (SNR)  and  to 
transfer  a  significant  portion  of  the  pump  wave  energy  to  an  incoherent  pedestal  remaining  after  compression. 
The  ratio  of  coherent  amplified  signal  energy  and  incoherent  superfluorescence  pedestal  energy,  i.e .,  signal-to- 
pedestal  ratio  (SPR),  has  been  measured  as  low  as  4:1  [2],  and  its  degradation  remains  a  limitation  to  pulse 
energy  scalability.  Due  to  the  intrinsic  difficulty  of  their  experimental  measurement  and  the  nonlinear  and 
quantum-mechanical  nature  of  the  physical  system,  the  saturation  dynamics  of  quantum  noise  in  OPCPA  have 
remained  controversial.  Their  understanding,  however,  is  requisite  for  optimization  of  OPCPA  performance. 

In  this  work,  we  have  conducted  a  numerical  analysis  of  saturation  dynamics  of  quantum-noise  contaminated 
OPCPA  employing  a  new  approach  which  captures  the  dynamics  of  energy  exchange  during  amplifier  saturation 
between  what  will  become  coherent  and  incoherent  components  of  the  electric  field  after  compression,  and  well 
reproduces  the  macroscopic  characteristics  observed  in  experiments  upon  comparison.  Our  approach  utilizes  the 
correspondence  between  the  evolution  equations  of  quantum  mechanical  expectation  values  and  the  classical 
Fokker-Planck  equations  even  for  nonlinear  systems,  when  the  nonlinearity  is  suitably  weak  [4].  This  allows  us 
to  accurately  compute  quantum-mechanical  expectation  values  from  stochastic  distributions  of  a  signal  and 
vacuum  noise  seeded  OPCPA  system  solved  numerically  from  the  classical  equations  of  motion. 
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Fig.  1  SNR  and  SPR  versus  crystal  length  for  a  2-pm  ultra-broadband  OPCPA  system.  Configurations  I,  II,  and  III 
correspond  to  a  system  with  signal  chirp  optimized  for  bandwidth,  efficiency,  or  efficiency-bandwidth  product, 
respectively,  as  defined  in  [5].  For  all  cases,  full  saturation  at  the  pump  pulse  peak  occurs  very  near  z  =  3  mm. 

We  find  several  distinguishing  features.  SNR  and  SPR  display  different  evolution  dynamics,  and  an  OPCPA 
has  well  defined  but  different  operating  points  for  maximum  suppression  of  noise-induced  fluctuations  or 
pedestal.  Beyond  these  operating  points  heavy  saturation  leads  to  large  excess  noise  that  can  be  enhanced  by 
orders-of-magnitude.  Comparison  of  this  analysis  to  a  2-jum  amplifier  system  [3]  yields  good  agreement.  In 
particular,  we  find  that  while  SNR  generally  improves  during  saturation  (i.e.,  shot- to-shot  fluctuations  are 
suppressed),  the  SPR  degrades,  especially  for  OPCPA  systems  with  a  high  degree  of  saturation  and  low  degree 
of  signal  chirp  chosen  to  optimize  bandwidth.  These  observations  provide  important  guidelines  for  the  design 
and  operation  of  low-noise  OPCPAs  and  insight  into  the  dynamics  of  quantum  noise  in  a  chirped  pulse  amplifier 
generally. 
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